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EXECUTIVE  SUMMARY 


The  Installation  Resjoration  program  (IRP)  underway  at  numerous  Air  Force 
bases>*ias  identified  several  sites  with  contaminated  soil  and  groundwater. 

This  subsurface  contamination  is  the  result  of  fuels,  cleaning  solvents,  and 
deqreasers  entering  the  subsurface  environment  from  accidental  spills,  leaking 
storage  tanks,  and  past  disposal  practices.  HQ  AFESC/RDVW  is  conducting 
research  aimed  at  developing  treatment  strategies  for  groundwater  cleanup,  and 
studying  the  fate  and  transport  of  contaminants  in  subsurface,  systems .  Many 
of  the  contaminants  of  concern  are  volatile  by  nature,  and  a  knowledge  of 
their  air-water  distribution  and  aqueous  solubility  is  needed  to  assess  the 
compounds  treatability  and  to  support  the  basic  laboratory  studies. 

The  objectives  of  this  research  were  to  develop  Henry's  law  constants  and 
aqueous  solubilities  as  a  function  of  temperature,  for  a  variety  of  organic 
compounds  of  Air  Force  concern  tTatrie  T).  Secondary  objectives  were  to 
determine  what  effect  mixed  organics,  in  an  aqueous  solution,  exhibit  on 
individual  Henry's  law  constants  and  evaluate  various  methods  used  to  predict 
Henry's  law  constants . 

This  report  documents  experimentally  determined  values  of  Henry's  law 
constants  and  aqueous  solubility  for  51  compounds  of  Air  Force  concern.  The 
report  is  presented  in  three  volumes.  Volume  I  contains  the  technical 
discussion  and  tabulated  values  of  Henry's  law  constants  and  aqueous 
solubilities.  Volume  II  and  III  contain  all  the  raw  data  and  the  fortran 
source  code  for  an  interactive  program  used  to  predict  the  chemical  parameters. 

Many  of  the  contaminants  of  concern  are  volatile  by  nature,  and  a 
knowledge  of  their  air-water  distribution  is  required  for  the  design  of 
treatment  processes  and  for  providing  insight  into  their  environmental  fate 
and  transport.  A  static  headspace  method  (Equilibrium  Partitioning  In  Closed 
Systems,  referred  to  as  EPICS)  was  used  to  measure  the  Henry's  law  constants, 
with  the  standard  batch  air  stripping  method  used  as  a  check. 

The  Henry's  constants  were  determined  as  a  function  of  tempertaure  from  10  to 
30  °C  (Table  11)  and  these  values  were  then  used  to  generate  temperature 
regression  equations  (Table  8).  Generally  speaking  the  EPICS'  results  from 
this  study  agree  well  with  other  published  results  (Table  T 2 ) . ^  However,  for 
many  of  the  compounds  reported  here,  confirmed  values  of  Henry's  constant  do 
not  exist  in  the  literature,  and  if  they  do,  values  are  rarely  reported  as  a 
function  of  temperature  with  rigorous  statistics. 

Solubility  data  for  organic  compounds  in  water  are  important  for 
environmental  studies  because  they  provide  fundamental  information  necessary 
to  predict  transport  in  aqueous  systems.  This  data  may  also  be  used  to 
predict  carbon  sorption  of  contaminants,  and  the  air-  or  stream-stri ppi nq 
behavior  for  a  given  compound.  The  aqueous  solubility  of  the  51  study 


compounds  were  determined  at  10,  20,  and  30°C  (Table  14).  Three  different 
methods  were  used,  but  the  majority  of  the  data  were  collected  using  a 
shake-flask  technique.  Although  the  solubilities  were  not  a  strong  function 
of  temperature  over  the  range  studied  (i.e.,  10-30  °C),  several  general 
trends  were  noted.  First,  the  solubility  of  the  halongenated  hydrocarbons 
increased  with  temperature.  Second,  the  solubility  of  the  substituted 
aromatic  hydrocarbons  increased  with  temperature.  Finally,  maxima  and  minima 
are  observed  for  a  wide  range  of  compounds  without  any  general  trend  that  can 
be  demonstrated  to  be  statistically  significant. 

Groundwater  contamination  is  often  characterized  by  the  presence  of 
several  different  contaminants  rather  than  one  single  compound.  For  this 
reason,  studies  were  conducted  to  determine  whether  the  presence  of  other 
compounds  would  affect  the  Henry's  law  constant  of  a  single  compound. 
Deviations  from  ideal  behavior  were  observed  (pg  52),  but  confirming 
experiments  were  not  performed.  Although  the  results  were  not  conclusive,  the 
project  team  believes  the  observed  interactions  were  real  and  reproducible. 

It  would  not  be  feasible  to  experimentally  determine  Henry's  law  constants 
for  all  chemical  compounds.  There  will  be  times  when  a  Henry's  law  constant 
is  needed  but  an  experimentally  determined  value  is  not  reported  and  the 
situation  does  not  permit  a  laboratory  study  to  determine  the  constant.  For 
this  reason,  a  technique  to  accurately  estimate  Henry's  constant  using  a 
minimum  of  physiochemical  properties  would  be  useful.  Three  different 
thermodynamic  techniques  for  correlating  experimental  Henry's  law  constants 

were  examined  (paqe  61).  The  techniques  were  examined  to  determine  their 
applicability  to  environmental  systems  and  their  predictive  capacity  for 
unmeasured  multicomponent  systems.  The  UNIFAC  method  proved  to  be  the  most 

effective  way  of  utilizing  tne  data  base  developed  durina  this  project.  A 
computer  algorithim  to  fit  the  current  data  to  a  new  environmental  UNIFAC 
binary  interaction  data  base  was  developed  and  a  portion  of  the  experimental 
data  collected  was  incorporated  into  this  new  data  base.  The  new  data  base 
creates  improvement  in  the  predictions  generated  by  UNIFAC  in  the  dilute 
concentration  regime  (Figures  13  through  16). 
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APPENDIX  C 


SUMMARY  OF  AQUEOUS  SOLUBILITY  MEASUREMENTS 


This  is  a  self-contained  document  with  its  own  internal  style, 
which  varies  from  our  format. 


1 


SOLUBILITY  MEASUREMENTS 


The  shake-flask  measurements  to  determine  solubility  required  several 
calibration  points  to  insure  their  statistical  validity.  In  this  study,  at 
least  five  independently  prepared  solutions  of  unknown  concentration  were 
calibrated  against  the  response  of  the  measurement  instrument.  i In  many,  if 
not  most,  previous  studies  serial  dilutions  of  one  sample  were  used  instead). 
In  the  case  of  aromatic  compounds  this  instrument  was  a  ultra  violet  (U.V. ) 
spectrometer.  The  calibration  samples  were  prepared  with  distilled  water  in 
gas  tight  amber  bottles  containing  no  headspace.  Since  water  was  employed  for 
U.V.  tests,  only  concentration  below  the  saturation  level  could  be  prepared. 

In  some  cases,  these  calibration  samples  might  be  too  far  below  the  measured 
U.V.  response  and  additional  samples  closer  to  the  true  solubility  limit  would 
be  prepared. 

The  data  was  logged  onto  a  Lotus ^  spreadsheet  and  a  linear  regression  of 
the  calibration  concentrations  versus  detector  response  determined.  The 
detector  response  for  the  saturated  solution  at  all  three  temperatures  couid 
then  be  entered,  and  the  corresponding  solubility  limits  determined  from 
extrapolation  of  the  linear  regression. 

In  the  case  of  the  gas  chromatograph  samples,  a  similar  scheme  was 
employed  using  integrator  peak  area  as  the  detector  response  instead  of  U.V. 
absorbance.  In  this  case,  however,  calibration  samples  were  prepared  in 
methanol,  so  calibration  concentrations  above  the  aqueous  solubility  limit 
might  be  prepared. 

Based  upon  the  correlation  coefficient  for  the  linear  regression  fit  on 
the  calibration  samples  (which  Is  shown  In  the  following  plots).  95  percent 
confidence  bands  can  be  placed  on  the  saturated  simple  concentrations.  These 
are  also  shown  on  the  following  plots.  Notice,  of  course,  that  the  better  the 
correlation  coefficient  the  narrower  the  bands.  Unfortunately,  in  most  cases 
(even  with  correlation  coefficients  in  excess  of  0.95).  little  statistical 
confidence  may  be  placed  in  the  temperature  dependency  of  the  solubilities 
measures.  As  a  consequence,  the  data  has  not  been  tabulated  explicitly  as  a 
function  of  temperature  dependence. 

The  nephelometry  results  are  not  listed  in  the  same  fashion  as  the  other 
results  due  to  differences  in  the  preparation  procedure.  However,  a 
discussion  of  the  statistics  of  nephelometry  is  contained  in  Section  IV  of  the 
main  text. 
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SOFTWARE  DOCUMENTATION 


This  Appendix  contains  the  Fortran  source  code  for  the  simplex  UNI FAC 
parameter  fitting  routine  and  the  interactive  program  for  calculating  Henry’s 
law  constants  and  aqueous  solubilities.  The  output  from  an  example  simplex 
fitting  run  is  included  for  the  interested  reader.  A  complete  Fortran  listing 
of  the  UN I FAC  VLE  and  LLE  binary  interaction  databases  is  also  presented  in 
this  Appendix  to  complete  the  software  documentation. 


EXAMPLE  PARAMETER  FIT  WITH  SIMPLEX  ALGORITHM 


******** ****** 
*  * 

*  CH2CQ/CCOH  * 

*  * 

************** 


ACETONE (11) -ETHANOL(20)  AT  305.15,  CAN. J. RES.  248,254,(1946) 
ETHANOL (20) -2, BUTANONE (12)  AT  1  ATM,  ZH.PRIKL.KHIM.  41,589,(1968) 
ISOPROPANOL (1) -2, PENTANONE( 16)  AT  1  ATM,  IEC  45,1803,(1953) 


COMPONENT  /  GROUPS 


1  2  17  15  22 

11  1  0  0  0  1 

20  0  0  1  0  0 

12  110  0  1 
16  1  2  0  0  1 

1  10  0  10 


GROUP  NO 

1 

2 

17 

15 

22 


GROUP  R 

.9011E+Q0 
.6744E+00 
.2106E+01 
. 1878E+01 
.1672E+01 


GROUP  Q 

.8480E+00 
.5400E+00 
.1972E+01 
. 1660E+01 
.  1488E+01 


GROUP  INTERACTION  PARAMETERS 


.OOOOE+OO 

•OOOOE+OO 

-.8793E+02 

-.8793E+02 

.2676E+02 


, OOOOE+OO 
.OOOOE+OO 
-.8793E+02 
- .8793E+02 
. 2676E+02 


.7375E+03 
.7375E+03 
.OOOOE+OO 
.OOOOE+OO 
. 1000E+01 


.7375E+03 
.7375E+03 
.OOOOE+OO 
.OOOOE+OO 
. 1000E+01 


.4764E+03 
.4764E+03 
. 1000E+01 
. 1000E+01 
.OOOOE+OO 


INITIAL  PARAMETERS  3  5 

100.000  100.000 
110.000  100.000 
100.000  110.000 
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SIMPLEX  ALGORITHM  RESULTS 


FINAL  PARAMETERS 
117.240  36.630 

FMIN=  .54424E-01  SD=  .70283E-G6 


TEMP  NUMBER  X  GAMEXP  GAMCAL  DEV 


305.15 

11 

.0750 

2.0516 

1.8340  - 

•10.6 

305.15 

20 

.9250 

1.0055 

1.0038 

-.2 

2 

305.15 

11 

.5000 

1.2187 

1.2088 

-.8 

305.15 

20 

.5000 

1.2304 

1.1903 

-3.3 

3 

305.15 

11 

.8000 

1.0356 

1.0330 

-.3 

305.15 

20 

.2000 

1.6686 

1.5984 

*4  ♦  2 

4 

305.15 

11 

1.0000. 

1.0000 

1.0000 

.0 

305.15 

20 

.0000 

2.3166 

2.1507 

-7.2 

5 

348.15 

20 

.2380 

1.4546 

1.5372 

5.7 

348.15 

12 

.7620 

1.0506 

1.0463 

-.4 

6 

347.55 

20 

.3970 

1.2497 

1.3013 

4.1 

347.55 

7 

347.55 

12 

.6030 

1.1154 

1.1306 

1.4 

20 

.6220 

1.0965 

1.1057 

.8 

347.55 

12 

.3780 

1.2825 

1.3382 

4.3 

8 

348.45 

20 

.7860 

1.0234 

1.0321 

.8 

348.45 

12 

.2140 

1.5091 

1.5754 

4.4 

9 

350.15 

20 

.9340 

1.0046 

1.0030 

-.2 

350.15 

12 

.0660 

1.7549 

1.8757 

6-9 

10 

371.37 

1 

.0755 

1.5221 

1.5656 

2.9 

371.37 

16 

.9245 

.9678 

1.0028 

3.6 

11 

361.71 

1 

.3720 

1.2395 

1.2481 

.7 

361.71 

16 

.6280 

1.0060 

1.0742 

6.8 

12 

358.59 

1 

.5845 

1.0974 

1.1070 

.9 

358.59 

16 

.4155 

1.1203 

1.2018 

7.3 

13 

356.76 

1 

.7815 

1.0257 

1.0298 

.4 

356.76 

16 

.2185 

1.2816 

1.4073 

9.8 
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DEBUG 


*****•*♦#■*  +  *#■*♦*#■*■*•**•**#****##****♦**♦***♦*♦*****#*********♦#♦***•*♦*** 


*  UNIFAC  PARAMETER  ESTIMATION  BY  MEANS  OF 

*  NELDER-MEAD'S  EXTENDED  SIMPLEX  MINIMIZATION  METHOD 


* 

* 


IMPLICIT  REAL* 8 (A-H,G-Z) 

DIMENSION  R< 10) ,  RS < 10),  QS < 10),  XL (10),  KG < 10 ) ,  NTEXT ( 40 ) ,  NUM ( 10 ) , 
&IPAR ( 4 ) ,  PAR ( 4 ) ,  DEV ( 100,  2 ) ,  JENS ( 10 ) , NU < 100, 2) 

DIMENSION  RI<76),  QI<76),  AI<40, 40), MAIN < 76) 

DIMENSION  X ( 5,  4 ) ,  F ( 5 ) ,  XB  <  4 ) ,  XS ( 4 ) ,  XM  <  4 ) , XE ( 4 ) , XX ( 4 ) , XR  <  4 ) , XK  <  4 ) 
COMMON  /U/  T ( 100 ) ,  NM ( 100,  2  > ,  XXX ( 100,  2 ) ,  GME ( 100,  2 ) ,  GMC ( 100,  2 ) , 
&GMR<100, 2), NNY<10, 10 ) , Q < 10 ) , A < 10, 10 ) 

DATA  NC,  NP/5,  6/ 

OPEN <6, FILE* 'SIMPLEX.  DAT',  ST ATUS= ' OLD ',  ACCESS* ' SEQUENTIAL ' , 
♦FORM* 'FORMATTED ' ) 

OPEN ( 7, FILE= 'SIMPLEX.  OUT ' ,  STATUS* ' NEW ' ,  ACCESS= ' SEQUENTIAL ' , 
♦FORM* ' FORMATTED ' ) 


READ  <  6,  IS )  NCOMP,  NG,  NOBS,  LAAF, NOIT, MDA 
16  FORMAT (2013) 


NCOMP*  NUMBER 
NG  *  NUMBER 
NOBS  *  NUMBER 
LAAF,  LAAF*1 

LAAF =2 

NOIT,  NQIT= 1 
N0IT=2 


MDA,  MDA*  1 
MDA  =  2 


OF  COMPONENTS 
OF  DIFFERENT  GROUPS 
OF  DATA  POINTS 

THE  EXPERIMENTAL  ACTIVITY  COEFFICIENTS 
ARE  READ 

THE  LOGARITHMS  TO  THE  EXPERIMENTAL 
ACTIVITY  COEFFICIENTS  ARE  READ 
PARAMETER  ESTIMATION  IS  PERFORMED 
NO  PARAMETER  ESTIMATION.  THE  ACTIVITY 
COEFFICIENTS  ARE  CALCULATED  BASED  ON 
THE  GIVEN  R,H  AND  A  MATRIX 
BUILT-IN  VLE  PARAMETERS  ARE  LOADED 
BUILT-IN  LLE  PARAMETER3  ARE  LOADED 


IF  <  NOIT. EQ. 2 )  GOTO  1301 
WRITE ( 7, 15) 

WRITE ( *, 15) 

15  FORMAT < 1H  ,'  PARAMETER  ESTIMATION') 

GOTO  1302 

1301  WRITE ( 7. 1303 ) 

WRITE ( * , 1303) 

1303  FORMAT < 1 H  CALCULATION  OF  ACTIVITY  COEFFICIENTS  BASED  ON  THE 

&GIVEN  Rr  Q,  AND  A  MATRIX') 

1302  CONTINUE 
DO  14  1  =  1,  16 
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READ (6, 13)  NTEXT 
WRITE <  *, 13  >  NTEXT 
14  WRITEC7, 13)  NTEXT 
13  FORMAT (40A2) 

READC6,  16)  <  NUM (  I ) ,  I  =  1, NCOMP ) 

NUM  GIVES  THE  NUMBERS  ATTACHED  TO 
THE  DIFFERENT  COMPONENTS 


DO  1304  1*1, NCOMP 
1304  READ  <6,  2)  (  NNY  (I,K),K  =  1,NG) 

2  FORMAT (4012) 

NNY ( I, K )  IS  THE  MATRIX  GIVING  THE  NUMBER  OF  GROUPS  OF 
KIND  K  IN  MOLECULE  I 

READ  <6, 16 )  <  KG (  K  ) ,  K  =  1,  NG  ) 

KG ( K )  IS  THE  NUMBER  ATTACHED  TO  GROUP  K 


DO  100  N=l, NOBS 

100  READ  (  6,  101 )  T(N),  <NU<N,  I),  XXX  <  N,  I ) ,  GME  (N,  I),  1*1,2) 

101  FORMAT ( F7. 2, 3 ( 14, 2F7. 4 ) ) 


N 

T  <  N ) 

NU<N, I)  * 

XXX  (  N,  I )  * 
GME  ( N,  I), 


DATA  POINT  NUMBER 
TEMPERATURE  IN  K 
NUMBER  ATTACHED  TO  COMPONENT  I 
<  NU<  N, I ) =NUM ( I ) ) 

LIQUID  MOLE  FRACTION 

LAAF*1  EXPERIMENTAL  ACTIVITY  COEFFICIENT 
LAAF=2  LOGARITHM  TO  THE  EXPERIMENTAL  ACTIVITY 
COEFFICIENT 


IF  <  LAAF-1 ) 620, 620, 821 

821  CONTINUE 

DO  822  N»  1, NOBS 
DO  622  1*1, 2 
GGG=GME<N, I) 

GME ( N, I ) =DEXP ( GGG ) 

822  CONTINUE 
820  CONTINUE 

IF  <  MDA. EQ. 1 )  CALL  UVLE ( RI, QI , AI , MAIN ) 
IF ( MDA. EQ. 2 )  CALL  ULLE < RI, QI, AI, MAIN ) 
DO  705  1=1, NG 
R  < I >  =RI <  KG  < I > ) 

Q ( I ) =QI <  KG  < I )  ) 

DO  705  IU=i, NG 

A (I,  IU)= AI<MAIN<KG< IU) ), MAIN <KG< I)  ) ) 
705  CONTINUE 

READ  <  6,  •  )  IOWN 

IF( IOWN. EQ. 0)  GOTO  707 
DO  704  IU=1, IOWN 
READ  <6,*  )  KOWNI 
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READ (6,  706)  R ( KOWNI ) , Q < KOWNI ) ,  < A (KOWNI, J), J*l, NG ) 

706  FORMAT ( 8F10. 4 ) 

704  CONTINUE 

707  CONTINUE 
C 

C  R<I)  IS  THE  GROUP  VOLUME  OF  GROUP  I 

C  Q<I>  IS  THE  GROUP  AREA  OF  GROUP  I 

C  Ad,  J)  IS  THE  GROUP  INTERACTION  PARAMETER 

C  BETWEEN  GROUPS  I  AND  J 

C  IOWN  IS  THE  NUMBER  OF  SETS  OF  USER -SUPPLIED 

C  R,  Q,  AND  A  <  I,  J )  DATA 

C  KOWNI  IS  THE  GROUP  ID  NUMBER  FOR  A  SET  OF 

C  USER-SUPPLIED  DATA 

C 

DO  7  1*1, NCOMP 
RS< I ) =0. DO 
QS(I)=Q. DO 
DO  8  J*l, NG 

RS< I ) =RS< I ) ♦  NNY < I,  J) *R( J) 

QS ( I ) =QS (  I ) +  NNY (I,  J) *Q ( J ) 

8  CONTINUE 

7  XL( I ) =5. DO*  <  RS( I ) -QS< I )  ) ~RS( I ) *1. DO 
WRITE <7, 5) 

WRITE  <  *,  5 ) 

5  FORMAT  < 1H  ,  '  COMPONENT  /  GROUPS',/) 

WRITE<7, 12)  (KG<I), 1*1, NG) 

WRITE( *, 12)  (KG<I), 1*1, NG) 

12  FORMAT  <17X,  1013) 

DO  10  1*1, NCOMP 

WRITE ( *, 6 )  NUM(I),  ( NNY ( I,  J ) ,  J=l, NG ) 

10  WRITE  <  7,  6  )  NUM(I),  <  NNY  (  I,  J  ) ,  J*  1,  NG  > 

6  FORMAT( 10X, 13, 4X, 1013) 

WRITE (7,3) 

WRITE ( *, 3) 

3  FORMAT ( //, '  GROUP  NO  GROUP  R  GROUP  Q', 

DO  11  1*1, NG 

WRITE ( *, 4 )  KG  <I),R(I),Q<I) 

11  WRITE ( 7, 4 )  KG(I),R(I),Q(I) 

4  FORMAT ( 4X, 13, 5X, 2E12. 4 ) 

WRITE (7, 1400) 

WRITE ( *, 1400) 

1400  FORMAT ( //,  '  GROUP  INTERACTION  PARAMETERS',/) 

DO  1401  1=1, NG 

WRITE< *, 1402)  ( A<I, J), J*l, NG) 

1401  WRITE ( 7, 1402  >  ( A < I, J ) , J  =  l, NG ) 

1402  FORMAT ( 10E12. 4 ) 

DO  102  N* 1, NOBS 
DO  102  1*1, 2 

DO  104  J=l, NCOMP 

IF(NU(N, I ) -NUM ( J ) >104, 103, 104 

103  NM ( N, I ) = J 
GOTO  102 

104  CONTINUE 
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102  CONTINUE 

CALL  PFAC3<RS,  QS, XL,  NOBS) 

IF ( NQIT-1 >830, £30, 831 
£31  CONTINUE 

CALL  PFAC4 <  NG, NOBS ) 

DO  £32  NR  =  1,N0BS 
DO  832  1=1, 2 

GMR ( NR,  I ) =  GMC <  NR,  I  > +GMR <  NR,  I ) 

GMR ( NR,  I ) =  DEXP (GMR  <  NR, I) ) 

£32  CONTINUE 
GOTO  833 
830  CONTINUE 

READ  <  6, 16 )  NPAR, KRIT, IDEN 
C 

C  NPAR  IS  THE  NUMBER  OF  PARAMETERS  TO 

C  BE  ESTIMATED 

C 

C  KRIT  DETERMINES  THE  OBJECTIVE  FUNCTION 

C  IF  KRIT=1,  THE  SUM  OF  THE  SQUARED  DIFFERENCES  BETWEEN 

C  THE  EXPERIMENTAL  AND  CALCULATED  ACTIVITY  COEFFICIENTS 

C  IS  MINIMIZED 

C  IF  KRIT=2,  THE  LOGARITHMS  OF  THE  ACTIVITY  COEFFICIENTS 

C  ARE  USED 

C 

C  IDEN  IS  THE  NUMBER  OF  IDENTICAL  PAIRS 

C  OF  INTERACTION  PARAMETERS 

C 

READ <6, IS)  < I PAR  <I), 1  =  1, NPAR ) 

IPAR  IS  THE  VECTOR  INDICATING  THE 
PARAMETERS  TO  BE  ESTIMATED  IN  THE 
A  MATRIX 

IF( IDEN) 950, 950, 951 
951  IDEN=2*IDEN 

READ<6,  IS)  <  JENS  <  J  ) ,  J  =  l, IDEN ) 

950  CONTINUE 

JENS  IS  THE  VECTOR  INDICATING  THE  IDENTICAL 
PAIRS  OF  PARAMETERS 

NN=NPAR+1 
N=NPAR 

READ ( 6, 400 )  ( X < 1, I ) , 1=1, NPAR ) 

XC1, I)=  ROW  OF  INITIAL  PARAMETERS 
NO  INITIAL  PARAMETER  MUST  BE  ZERO 


SA=1. D-6 


SA  IS  THE  STANDARD  ERROR  AS  DEFINED  BY 
NELDER-MEAD 
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DO  201  J=2, MN 

DO  201  1  =  1, N 

IF< J-I -1)202, 203, 202 

203  X< J,  I)=l. 1D0*X(1, I) 

GOTO  201 

202  X<  J,  I)=X! 1,  I  ) 

201  CONTINUE 

WRITE  !  7,  300  )  <  IPAR  <I),I  =  1,N> 

WRITE! *, 300)  < IPAR (I), 1=1, N) 

300  FORMAT!//,'  INITIAL  PARAMETERS 413/ ) 

DO  204  J=l, NN 

WRITE  C  * , 400 )  <X< J, I), 1  =  1, N) 

204  WRITE! 7, 400)  ! X < J, I ) , I =1, N ) 

DO  1  J  =  l,  NN 

DO  21  1  =  1,  N 
21  XX< I >  =X<  J, I) 

CALL  FMIN  (  NPAR,  IPAR,  PAR,  NOBS,  NG,  XX,  FF,  KRIT,’ JENS,  IDEN  ) 
1  F<  J ) =FF 
NF  =  NN 

NF  IS  THE  NUMBER  OF  CALCULATIONS  OF  F 

ALFA=i. DO 
BETA=0. 5D0 
GAMMA=2. DO 
ITER=0 
JPR*0 

400  FORMAT C8F10. 3) 

ESTIMATION  OF  THE  LOWEST  VALUE  OF  F=FB 

25  FB=F ( 1 ) 

DC  98  1=1, N 
98  XB ( I ) =X  < 1,  I) 

JB  =  1 

DO  31  J =2, NN 
IF ( FB-F ( J ) >31, 31, 108 
108  FB=F  <  J ) 

JB= J 

DO  41  1=1, N 
41  XB  <  I )  =X  !  J,  I) 

31  CONTINUE 

ESTIMATION  OF  THE  HIGHEST  VALUE  OF  F=FS 

FS=F ! 1 ) 

DO  51  1=1, N 
51  XS< I ) =X  < 1, I > 

JS  =  1 

DO  61  J  =  2,  NN 
IF ( FS-F ! J ) ) 111, 61, 61 
111  FS  =  F  <  J ) 

JS  =  J 
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DO  7 1  1*1, N 
71  XS<I>=X<J,I) 

61  CONTINUE 

CALCULATION  OF  THE  CENTROID  XMCI)  QF  POINTS 
EXCLUDING  XS(I) 

DO  81  1*1, N 
81  XM<I)*-XSa> 

DO  9  J=i, NN 
DO  122  1  =  1,  N 
122  XM  < I ) *XM  < I > ♦X  <  J,  I> 

9  CONTINUE 
DO  121  1  =  1,  N 
121  XM  < I ) =XM  < I > /DBLE ( N ) 

REFLECTION 

DO  131  1  =  1,  N 

131  XR ( I ) =XM ( I ) ♦ALFA* ( XM ( I ) -XS< I ) > 

CALL  FMIN  <  NPAR,  IPAR,  PAR,  NOBS,  NG,  XR,  FR,  KRIT,  JENS, IDEN ) 
NF=NF*1 

EXPANSION 

IF ( FR-FB ) 141, 151, 151 
141  DO  161  1=1, N 

161  XE  ( I )  =XM  ( I )  ♦GAMMA  *  ( XR  C I )  -XM  ( I )  ) 

CALL  FMIN  <  NPAR,  IPAR,  PAR,  NOBS,  NG,  XE,  FE,  KRIT,  JENS,  IDEN  ) 
NF=NF+1 

IF(FE-FB) 17, 18, 18 
17  DO  19  1=1, N 
X(  JS,  I ) =XE  < I ) 

19  XS ( I ) =XE  < I ) 

F( JS) =FE 

CALCULATION  OF  THE  HALTING  CRITERION 

27  FM=0. DO 

DO  20  J=l, NN 

20  FM  =  FM  +  F  <  J ) 

FM=FM/DBLE ( NN ) 

FRMS=0. DO 

DO  22  J*l, NN 

22  FRMS* (  F  ( J  ) -FM ) *  *2  +  FRMS 
RMS=DSQRT  <  FRMS/DBLE ( N )  ) 

ITER= ITER+1 
JPR= JPR+1 

IF ( ITER-200 ) 500, 500, 23 
500  CONTINUE 

IF<  JPR-D902,  902,  903 
903  CONTINUE 

IF  <  JPR-6 ) 9Q1, 904, 904 
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904  JPR  =  1 
902  CONTINUE 

WRITE  <  7, 107) 
WRITE ( *, 107) 
107  FORMAT  <  / /,  ' 
4E ' ,  150 
WRITE  <  7, 109) 
WRITE <  *,  109) 
109  FORMAT ( ' 

WRITE (7,  400) 
WRITE <  *, 400) 
WRITE  (  7, 106) 
WRITE  <  *,  106) 
106  FORMAT <1H  ,  ' 
901  CONTINUE 

IF ( RMS-SA ) 23 


ITER, NF 
ITER, NF 
ITERATION',  14, 


PARAMETERS' ) 

<X<  JS,  I),  1*1,  N  ) 
(X(  JS,  I),  1  =  1,  N) 
F< JS), RMS 
F  <  JS ) , RMS 
FMIN*  ' ,  E14.  5,  ' 

23,  25 


NUMBER  OF  CALLS  FOR  THE  SUBROUTIN 


SD= ' , E14. 5) 


NEW  SIMPLEX 
FE  GREATER  THAN  FB 


18  DO  26  1=1, N 
XCJS,  I ) =XR  < I ) 

26  XS ( I ) =XR ( I ) 

F ( JS ) =FR 
FS=FR 
GOTO  27 

NEW  SIMPLEX 
Fp  GREATER  THAN  FB 

111  DO  30  J= 1, NN 

IF ( J- JS ) 28, 30, 28 
28  IF(FR-F< J) ) 18, 18, 30 
30  CONTINUE 

IF ( FR-FS ) 91, 91, 32 
91  DO  33  1  =  1,  N 
X<  JS,  I )  *XR  <  I ) 

33  XS ( I ) =XR ( I ) 

F  <  JS ) =FR 
FS*FR 

32  DO  34  1  =  1,  N 

34  XK  < I ) =XM  < I ) +BETA  * ( XS ( I ) -XM  < I ) ) 

CALL  FMIN <  NPAR,  IPAR,  PAR,  NOBS,  NG,  XX, FK, KRIT, JENS, IDEN ) 
NF=NF-1 


NEW  SIMPLEX 
AFTER  CONTRACTION 


IF <  FK-FS ) 35, 35, 36 
35  DO  37  1=1, N 
X(  JS,  I  )  =XK  <  I ) 

37  XS ( I ) =XK ( I ) 

F( JS) =FK 
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FS=FK 
GOTO  27 

36  DO  3S  J=l, HH 
DO  39  1=1, N 

39  X(  J,  I)  =  (X<  J,  I )  +  XB<  I )  >/2.  DO 
38  CONTINUE 
GOTO  27 

23  WRITE ( *, 905) 

WRITE  <  7, 905) 

905  FORMAT (  //,  '  FINAL  PARAMETERS') 

WRITE  (  7,  400  >  <X<  JS,  I),  1  =  1, N> 

WRITE  (  *,  400  >  <X<  JS,  I),  1  =  1, N) 

WRITE<7, 106)  F ( JS ) , RMS 

WRITE<  *, 106 >  F ( JS ) , RMS 
833  CONTINUE 

DO  906  N=i, NOBS 
DO  906  1  =  1,  2 

DEV ( N,  I )  =  ( GMR ( N,  I)-GME<N,  I) )*1Q0. DO/GME  <  N, I) 

906  CONTINUE 
WRITE<7, 207) 

WRITE (  *,  207) 

207  FORMAT  (.//,’  TEMP  NUMBER  X  GAMEXP  GAMCAL  DEV',/) 

DO  222  N=i, NOBS 
WRITE<7,  16)  N 
WRITEC*,  16)  N 
DO  223  1  =  1,2 

WRITE  (  *,  96)  T(N),  NU(N,  I),  XXX(N,  I  >  ,  GME  (  N,  I ) ,  GMR  <  N,  I ) ,  DEV  ( N,  I  ) 

223  WRITEt  7,  96)  T  <  N ) ,  NU  <  N,  I ) ,  XXX  ( N,  I ) ,  GME<  N,  I ) ,  GMR  ( N,  I),DEV(N,  I) 

96  FORMAT ( F8.  2,  15,  FS.  4,  2F14.  4,  F6.  1 ) 

222  CONTINUE' 

IF ( NOIT- 1 ) 835,  835,  836 

835  CONTINUE 
WRITE ( 7, 1400) 

WRITE< *,  1400) 

DO  112  1=1, NG 

WRITEC  *,  1402)  C  A  <  I,  J),  J  =  i,  NG  > 

112  WRITE(7, 1402)  ( A ( I, J ) , J=i, NG ) 

WRITE (7,  996) 

WRITE < *,  996) 

996  FORMAT <//, '  THE  SUM  OF  THE  SQUARED  DIFFERENCES  BETWEEN  THE*) 

IF ( KRIT- 1 ) 997,  997,  998 

998  WRITE <  *, 995 ) 

WRITE ( 7, 995) 

995  FORMAT ('  LOGARITHMS  TO  THE') 

997  WRITE( *, 999 ) 

WRITE ( 7,  999) 

999  FORMAT < '  EXPERIMENTAL  AND  CALCULATED  ACTIVITY  COEFFICIENTS  IS 

£  MINIMIZED' > 

836  CONTINUE 
STOP 
END 
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•  •♦•♦♦"ft****-***#**-**********##**********##**************************** 


SUBROUTINE  PFAC3 


* 


SUBROUTINE  PFAC3  CALCULATES  THE  COMBINATORIAL  PART 
OF  THE  ACTIVITY  COEFFICIENTS 


SUBROUTINE  PFAC3 ( RS, QS, XL, NOBS  > 


IMPLICIT  REAL* 8 ( A-H,  O-Z ) 

DIMENSION  THETA ( 2  > ,  PHI <  2 ) ,  RS  < 10 ) , QS  < 10  > 

COMMON  /U/  TdOO),  NMCIOO,  2),  XXXdOO,  2>, 

&GMR  (100,  2),  NNY(10,  10  ) ,  Q  <10  ) ,  A  < 10,  1.0  > 

DO  3  N= 1, NOBS 
SQ-O. DO 
SR=0. DO 
SXL=Q. DO 
DO  2  1=1, 2 
J=NM  <  N,  I) 

SXL=SXL+XL ( J ) *XXX <  N,  I ) 

SQ*SQ*QS  <  J ) *XXX <  N,  I ) 

2  SR=SR*RS< J) *XXX(N,  I) 

DO  3  1  =  1,2 
J«NH<N, I) 

THETA  <  I )  =QS  (  J  )  /SO 
PHI < I ) =RS<  J ) /SR 

GMC<N,  I )  =DLOG  <  PHI  (  I )  )  +-5.  DO*QS  ( J  )  *DLOG  ( THETA  <  I )  /PHI  <  I )  )*XL(  J) 
&-PHI < I ) *SXL 

3  CONTINUE 
RETURN 
END 


, XL (10) 

GME  ( 100,  2 ) ,  GMC  (100,  2 ) , 


» 

*  .  SUBROUTINE  PFAC4 


*  SUBROUTINE  PFAC4  CALCULATES  THE  RESIDUAL  PART 

*  OF  THE  ACTIVITY  COEFFICIENTS 

* 


SUBROUTINE  PFAC4 ( NG, NOBS ) 


IMPLICIT  REAL*a< A-H, O-Z) 

DIMENSION  GMOLdO) ,  ATET(IO),  ANYK  ( 10  ) ,  BNYK  <10  ) ,  GK  (10,  3),  PdO,  10) 
COMMON  /U/  TdOO),  NMdOO,  2),  XXXdOO,  2),  GMEdOO,  2),  GMCdOO,  2), 
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&GMR( 100,  2  > ,  NNY ( 10,  10 ) ,  Q < 10 > , A < 10, 10 > 

CALCULATION  OF  THE  PSI  MATRIX 

DO  250  NR*1, NOBS 
DO  7  -1  =  1,  NO 
DO  7  J=l,  NG 

7  PfI,  J)=DEXP(-A< I,  J)/T(NR) ) 

CALCULATION  OF  GROUP  MOLE  FRACTIONS 

DO  105  11  =  1,3 

IF ( II-2 ) 100,  100, 101 

100  SNYK=0. DO 

PURE  COMPONENT 

J=NM  <  NR, II) 

DO  12  K=l, NG 

12  SNYK=SNYK>DBLE<NNY< J, K5 > 

DO  13  K  =  l, NG 

13  GMOL  <  K  >  =DBLE  <  NNY (  J,  K  >  ) /SNYK 
GOTO  102 

101  SNYK=Q. DO 

MIXTURE 

DO  2  1*1,2 
J=NM(NR, I) 

DO  2  K  =  1,  NG 

2  SNYK=SNYK+DBLE ( NNY< J,  K  > >*XXX(NR,  I) 

DO  3  K*l,  NG 

GNYK=Q.  DO 
DO  4  1  =  1,  2 
J=NM ( NR, I ) 

4  GNYK  =  GNYK*DBLE ( NNY  <  J, K )  >*XXX<NR,  I) 

3  GMOL ( K ) =GNYK/SNYK 

CALCULATION  OF  GROUP  AREA  FRACTIONS 

102  SNYK=Q. DO 
DO  5  K  =  l,  NG 

5  SNYK=SNYK*Q  <  K ) *GMOL  <  K ) 

DO  6  K=i,  NG 

6  ATET  t  K ) =0  <  K ) *GMQL ( K ) /SNYK 

CALCULATION  OF  GAMMA  K 


DO  9  K  =  i,  NG 
ANYK(K) =Q. DO 
DO  10  M=l, NG 
SNYK=0. DO 
DO  8  N*l,  NG 
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S  SNYK=SNYK*ATET(N) *P(N, «) 

10  ANYK (K)=ATET(M)*P(K, W > /SNYK+ANYK <K ) 

BNYK ( K ) =0. DO 

DO  11  M  =  1,NG 

11  BNYK(K)  =BNYK(K)  *ATET( M>  *P<  M,  K) 

BNYK  C.K  >  =DLQG  <  BNYK  <  K  >  ) 

9  GK ( K,  II)=Q(K)*(1. DO-BNYK(K) -ANYK(K) ) 

105  CONTINUE 

DO  201  1=1, 2 
J=NM(NR, I) 

SNYK=0. DO 
DO  200  K=l, NG 

200  SNYK=SNYK*DBLE  C  NNY <  J,  K ) ) *  <  GK  <  K,  3 ) -GK ( K,  I > ) 

201  GMR ( NR, I >  =SNYK 
250  CONTINUE 

RETURN 

END 


SUBROUTINE  FMIN 

SUBROUTINE  FMIN  CALCULATES  F  (THE  OBJECTIVE  FUNCTION) 
AS  A  FUNCTION  OF  A  SET  OF  PARAMETERS 


SUBROUTINE  FMIN ( NPAR,  IPAR,  PAR,  NOBS,  NG, XX,  FF,  KRIT, JENS,  IDEN ) 


IMPLICIT  REAL*B< A-H, 0-2) 

DIMENSION  IPAR (4), PAR (4), XX (4), JENS (10) 

COMMON  /U/  T(100),  NMC100,  2),  XXX ( 100,  2),  GME(100,  2),  GMC(100,  2), 
&GMR ( 100,  2),  NNY (10, 10 ) , Q < 10 ) , A < 10, 10 ) 

DO  2  1  =  1,  NPAR,  2 
KI  =  IPAR ( I ) 

KJ  =  IPAR( 1*1 ) 

A(KI,KJ)  =  XX  ( I ) 

2  A ( KJ, KI ) =XX (1*1) 

IF ( IDEN ) 9, 9, 8 

8  KKI = JENS ( 1 > 

KKJ= JENS( 2) 

DO  7  J=3, IDEN, 2 
IKI = JENS ( J ) 

IKJ= JENS ( J*1 ) 

A  < IKI,  IK J ) =A<  KKI, KKJ ) 

A ( IK J,  IKI ) = A ( KKJ, KKI) 

7  CONTINUE 

9  CONTINUE 
CALL  PFAC4 ( NG, NOBS ) 

DO  200  NR=1, NOBS 
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DO  200  1  =  1,2 

GMR  <  NR,  1 >  =GMC ( NR, I ) +  GMR ( HR,  I > 

GMR  (  NR,  I )  =  DEXP  ( G21R  (  NR,  I  )  ) 

200  CONTINUE 
FF =0. DO 
DO  3  .  N= 1, NOBS 
DO  3  1  =  1,2 
IF  <  KRIT-1 ) 10, 10, 20 
10  FF»FF«-<GMR<N, I>-GME(N, I> > 

GOTO  3 

20  GCAL=G«R(N, I) 

G£XP=GMECN, I) 

FF-FF+ < DLQG ( GCAL ) -DLQG ( GEXP >  >  *  *2 
3  CONTINUE 
RETURN 
END 


* 


UVLE  CONTAINS  BUILT-IN  UNIFAC  VLE -PARAMETERS 


SUBROUTINE  UVLE ( RI,  QI, AI, WAIN ) 
C  . . . * 

c 

C  THE  MAIN  GROUPS  ARS: 


c 

1 

CH2  .  . 

2 

.  3 

ACH . 

.  4 

ACCH2  . . . 

.  5 

QH . 

c 

6 

CH30H 

7 

H20 . 

.  a 

ACOH  . . . . 

.  •  9 

CH2C0  . . . 

10 

CHO . 

c 

11 

ccao  . 

•  •  •  • 

12 

HCOO  .... 

13 

CH20  . . . . 

14 

CNH2  .... 

15 

CNH . 

c 

16 

<C)3N 

•  •  «  • 

17 

ACNH2  . . . 

IS 

PYRIDINE 

19 

CCN . 

20 

COQH . 

c 

21 

CCL  .  . 

•  •  •  • 

22 

CCL2  .... 

23 

CCL. 3  •  •  •  • 

24 

CCL4  .... 

25 

ACCL . 

c 

26 

CN02  . 

•  •  «  • 

27 

ACN02  . . . 

28 

CS2 . 

29 

CH3SH  . . . 

30 

FURFURAL  . 

c 

31 

DOH  .  . 

32 

I . 

33 

BR . 

34 

C--C  .... 

35 

DMSO  ...  . 

c 

c 

36 

ACRY  . 

37 

CLCC  • • • • 

38 

ACF . 

39 

DMF . 

40 

CF2 . 

c 

c 

THE  SUB  GROUPS  ARE: 

c 

1 

CH3  .  . 

2 

CH2 . 

.  3 

CH . 

.  4 

c. . 

.  5 

CH2=CH 

c 

6 

CH  =  CH 

•  *  t  t  t 

7 

CH2=C  . . . 

.  a 

CH=C  . . . . 

.  9 

c=c . 

10 

ACH . 

c 

11 

AC  *  •  • 

•  •  •  * 

12 

ACCH3  . . . 

13 

ACCH2  . . . 

14 

ACCH  .... 

15 

OH . 

c 

16 

CH30H 

•  •  •  * 

17 

H2Q . 

18 

ACQH  .... 

19 

CH3C0  . . . 

20 

CH2CO  . . . . 

c 

21 

CHO  .  . 

22 

CH3C00  . . 

23 

CH2C00  . . 

24 

HCOO  .... 

25 

CH30 . 

c 

26 

CH20  . 

t  •  *  • 

27 

ch-q 

28 

FCH20 

29 

CH3NH2  . . 

30 

ruowu*5 

c 

31 

CHNH2 

•  •  »  • 

32 

CH3NH  . . . 

33 

CH2NH  . . . 

34 

CHNH  .... 

35 

CH3N . 

c 

36 

CH2N  . 

•  •tt 

37 

ACNH2  . . . 

28 

C5H5N  . . . 

39 

C5H4N  . . . 

40 

C5H3N  . . . . 

c 

41 

CH3CN 

•  •  •  * 

42 

CH2CN  . . . 

43 

COQH  .... 

44 

HCOOH  . . . 

45 

CH2CL  . . . . 

c 

46 

CHCL  . 

•  *  *  * 

47 

CCL . 

48 

CH2CL2  . . 

49 

CHCL2  . . . 

50 

CCL2  . . . . , 

c 

51 

CHCL3 

52 

CCL3 _ 

53 

CCL4  .... 

54 

ACCL  .... 

55 

CH3N02  . . . 

60 


c 

56 

CH2N02  . . . 

57 

CHN02  . . . 

58 

ACN02  . . . 

59 

CS2  . 

.  .  .  .  60 

CH3SH  . . . . 

c 

61 

CH2SH  .... 

62 

FURFURAL 

S3 

( CH2QH ) 2 

64 

I  ... 

.  .  .  .  65 

BR . 

c 

66 

CH  =  ~C _ 

67 

C--C  .... 

68 

DflSQ  .... 

69 

ACRY 

.  .  .  .  70 

CL  <  C  =  ~C )  . 

c 

71 

ACF  ...... 

72 

DMF-i  . . . 

73 

DMF-2  . . . 

74 

CF3  . 

.  .  .  .  75 

CF2 . 

C  76  CF. . 

C 

. . . . . ***************.**************************« 

c 

IMPLICIT  REAL *8  <  A-H, 0-2 ) 

DIMENSION  AX <  40,  40 ) ,  RR < 76 )  ,  QG  < 76 ) ,  RI ( 76 ) ,  01 ( 76 ) ,  MAINSG ( 76 ) , 

*  MAIN < 76) 

DIMENSION  Ai <40),  A2<40>,  A3 (40),  A4<40) ,  A5<  40), A6 ( 40) , A7 < 40 ) , AS<40', 

*  A  9 <  40 > , A10(40), Ali<40),  A12<4G>,  A13(4Q>,  A14<40>,  A15 ( 40 ) , A16 < 40 ) , 

*  A17  <  40 ) ,  A1S<40), A19<40>,  A20<40>,  A21 ( 40 ) ,  A22 ( 40 > ,  A23 \ 40 ) ,  A24 ( 40 ) , 
*A25( 40 ) ,  A26<40),  A27C4Q),  A2S<40>,  A29 < 40 ) ,  A30 ( 40 ) ,  A31 ( 40 ) , A32 ( 40 ) , 

* A33 ( 40 ) ,  A34 <  40 ) ,  A35 ( 40 ) ,  A36 ( 40 ) ,  A37 ( 40 ) ,  A38  <  40 ) ,  A39  <  40 ) ,  A4C ( 40  > 

DATA  MAINSG/  4*1,  5*2,  2*3,  3*4,  1*5,  1*6,  1*7,  1*8,  2*9,1*10, 

*  2*11, 1*12, 4*13, 3*14,  3*15,  2*16,  1*17,  3*18,  2*19,  2*20, 

*  3*21,  3*22,  2*23,  1*24,  1*25,  3*26,  1*27, 1*28, 2*29,  1*30, 

*  1*31, 1*32,  1*33,  2*34,  1*35,  1*36,  1*37,  1*38,  2*39, 3*40/ 

DATA  RR/O.  9011,  0.  6744,  0,  4469,  O.  2195,  1.  3454,  1.  1167,  1.  1173,  0.  8886, 

*0.  6605,  O.  5313,  0.  3652,  1.  2663,  1.  0396,  0.  8121,  1.  O,  1.  4311,  0.  92,  0.  8952, 
*1.  6724,  1.  4457,  0.  9980,  1,  9031,  1. 6764,  1.  2420,  1.  1450,  0.  9183,  0.  6908, 

*0.  9183,  1.  5959,  1.  3692,  1.  1417,  1.  4337,  1.  2070,  0.  9795,  1.  1865,  0.  9597, 

*1.  0600,  2.  9993,  2.  8332,  2.  6670,  1.  8701,  1.  6434,  1.  3013,  1.  5280,  1.  4654, 

*1.  2380,  1.  0060,  2.  2564,  2.  0606,  1. 8016,  2.  8700,  2.  6401,  3.  3900,  1.  1562, 

*2.  0086,  1. 7818,  1.  5544,  1. 4199,  2.  0570,  1.  8770,  1.  6510,  3.  1680,  2.  4088, 

*1.  2640,  O.  9492,  1.  2920,  1.  0613,  2.  8266,  2.  3144,  0.  7910,  O.  6948,  3.  0856, 

*2.6322,1.4060,1.0105,0,6150/ 

DATA  QQ/O.  848,  0.  540,  0.  228,  O.  000,  1.  176,  0.  867,  O.  988,  0.  676,  0.  485, 


►0.  400,  0.  120 

,  0.  968,  0.  660,  0.  348, 

1.  200,  1.  432 

,  1.  400,  0. 

680, 

1. 488, 

►1.  180,  0.  948 

,  1.  728,  1. 420,  1.  188, 

1.  088,  0.  780 

,  0.  468,  1. 

100, 

1. 544, 

►1.  236,  0.  924 

,  1. 244,  0.  936,  0.  624, 

0.  940,  0.  632 

,  0.  816,  2. 

113, 

1.833, 

►1.  553,  1.  724 

,  1.  416,  1. 224,  1.  532, 

1.  264,  0.  952 

,  0.  724,  1. 

988, 

1.684, 

>1.  448,  2.  410 

,  2.  184,  2.  910,  0.  844, 

1.  868,  1.  560 

,  1.  248,  1. 

104, 

1.  650, 

►  1.  676,  1.  368 

,  2.  481,  2.  248,  0.  992, 

0.  832,  1.  088 

,  0.  784,  2. 

472, 

2.  052, 

*0.  724,  0.  524 

,  2.  736,  2.  120,  1.  380, 

0.  920,  0.  460/ 

DATA  Al/ 

0.  000, 

-200.  000, 

61.  130, 

76. 500, 

986. 500, 

►  697.  200, 

1318. 000, 

1333. 000, 

476. 400, 

677.  000, 

►  232. 100, 

741. 400, 

251. 500, 

391. 500, 

255.  700, 

►  206. 600, 

1245. 000, 

287. 770, 

597. 000, 

663. 500, 

►  35.  930, 

53. 760, 

24. 900, 

104.  300, 

321. 500, 

►  661.500, 

543.  000, 

153. SCO, 

184.  400, 

354.  550, 

►  3025. 000, 

335. 800, 

479.  500, 

298. 900, 

526. 500, 

►  689. OOC, 

-0.  505, 

125.  800, 

485.  30C, 

-2.  859/ 

DATA  A2/ 

2520. 000, 

0.  000, 

340.  700, 

4102.  000, 

693. 900, 

-  1509.000, 

634. 000, 

547.  40G, 

524. 500, 

9000.  000, 

/  ,*►  »  — *  w  0 

468. 700, 

5<2Q  TOO 

w  w  ■  w  W  W  p 

ioc  nnr> 

«•  '  w  *  W  W  W  f 

►  658.  800, 

9000. 000, 

9000. 000, 

405. 900, 

730. 400, 

►  99.  810, 

337. 100, 

4584. COO, 

5831. 000, 

959.  700, 

►  542.  100, 

9000. 000, 

76.  302, 

9000. 000, 

9000.  000, 

►  9000. 000, 

9000. 000, 

9000. 000, 

523. 600, 

9000.  000, 

►  9000. 000, 

237. 30C, 

9000. 000, 

320. 400, 

9000.  000/ 
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DATA  A3/ 

-11.  120, 

-94. 780, 

♦ 

637. 350, 

903.  800, 

1229. OOO, 

♦ 

5.  994, 

9000. OOO, 

32.  140, 

♦ 

90. 490, 

668.  200, 

-4.  449, 

* 

-18.  810, 

-144.  400, 

-231.  900, 

# 

168. OOO, 

194.  900, 

52.  070, 

* 

210. 400, 

113.  300, 

-13.  590, 

* 

9000. OOO, 

69.  110, 

389.  300, 

DATA  A4/ 

-69. 700, 

-269.  700, 

* 

603. 250, 

5695. OOO, 

884. 900, 

* 

5688. OOO, 

9000. OOO, 

213.  100, 

* 

23. 500, 

764.  700, 

52. 800, 

♦ 

-114. 100, 

9000. 000, 

-12. 140, 

* 

3629. OOO, 

4448. OOO, 

-9.  451, 

* 

4975.  OOO, 

9000. OOO, 

-171.  300, 

♦ 

9000- OOO, 

9000. OOO, 

101.  400, 

DATA  AS/ 

156. 400, 

8694. OOO, 

* 

-137, 100, 

353. 500, 

-259.  700, 

♦ 

101. OGO, 

193. 100, 

28. 060, 

* 

-323. OOO, 

-348. 200, 

170.  OOO, 

* 

75.  620, 

-112. 100, 

-98.  120, 

* 

61. 110, 

157. ICO, 

477. OOO, 

♦ 

-318. 900, 

313.  500, 

133.  400, 

* 

9000. OOO, 

253. 900, 

44.  780, 

DATA  A6/ 

16. 510, 

-52. 390, 

» 

0.000, 

-181 . OOO, 

-101.  700, 

* 

-1G. 720, 

193. 400, 

-180. 600, 

♦ 

53.  900, 

335. 500, 

580. 500, 

» 

-38. 320, 

-102. 500, 

-139.  400, 

* 

75.  140, 

9000. OOO, 

-31. 090, 

# 

9000. OOO, 

9000. OOO, 

9000.  OOO, 

• 

9000. OOO, 

-21. 220, 

-48.  250, 

DATA  A7/ 

300. QOQ, 

692. 700, 

289. SOO, 

0.  OOO, 

324. 500, 

14.  420, 

9000. OOO, 

540. 500, 

304. OOO, 

213.  OOO, 

459.  OOO, 

325. 400, 

3:  .  400, 

353. 700, 

220. 600, 

399. 500, 

887. 100, 

0.  OOO, 

9000. OOO, 

9000. OOO, 

160. 800, 

9000. 000, 

9000. OOO, 

DATA  AS/ 

275. 800, 

1665.  OOO, 

-265. 200, 

-601. 800, 

0.  OOO, 

-449. 400, 

9000. OOO, 

9000.  OOO, 

9000. OOO, 

9000.  OOO, 

-305. 500, 

9000. OOO, 

9000.  OOO, 

9000.  OOO, 

9000. 000, 

9000. COO, 

9000. OOO, 

-687. 100, 

9000.  OOO, 

9000. OOO, 

annn  nnn 

w wW  W  1  W  W W  f 

9000.  OOO, 

anpn  nnn 

DATA  A9/ 

26.  760, 

-82.  920, 

* 

108. 700, 

472.  500,, 

-133.  100, 

♦ 

-213. 700, 

9000. 000, 

5.  202, 

# 

9000. OOO, 

937. 900, 

165.  100, 

• 

-191. 700, 

-284. OOO, 

-354. 600, 

0.000,  167.000,  636.100, 

25. 770,  9000. 000, 

161.700,  122.  S00, 

212.  500,  537.  400, 

3. 000,  538. 200, 

-10.430,  -64.690, 

9000.000,  169.900, 

245.  600,  9000.  000/ 

-146.800,  0.000,  803.200, 

-52. 100,  9000. 000, 

9000.000,  -49.290, 

6096.  000,  603.  800, 

-141.300,  -126.900, 

9000. 000,  -20. 360, 

9000. 000,  4284. 000, 

5629. 000,  9000. 000/ 

89.600,  25.820,  0.000, 

84.000,  441.800, 

83. 020,  42. 700, 

6.712,  199.000, 

143.100,  287.800, 

147.500,  -120.500, 

9000. 000,  -202. 100, 

-143.900,  9000.000/ 

-50.000,  -44.500,  249.100, 

23. 390,  306.  400, 

359. 300,  266. 000, 

36. 230,  -289. 500, 

-67.800,  17.120, 

37. 840,  9000. 000, 

9000. 000,  -399. 300, 

-172.400,  9000.  OGO/ 

362.300,  377.600,  -229.100, 

-195.400,  -257.300, 

48.890,  168.000, 

112.600,  -14.090, 

497.500,  678.200, 

9000. ooo,  iaa. ooo, 

9000.000,  -139.000, 

319.000,  9000.000/ 

25.340,  244.200,  -451.  60C, 

-356.100,  9000.000, 

9000. OOO,  9000. OOO, 

9000. OOO,  9000. OOO, 

1827.  COO,  9000.  OOO, 

9000. OOO,  9000. OOO, 

9000.000,  9000.000, 

9000.  OOO,  9000.000/ 

140.100,  365.800,  164.500, 

0.000,  -37.360, 

9000. COO,  9000. OOO, 

481.700,  669.400, 

-39.200,  174.500, 


*  137.500,  9000.000,  216.100, 

*  9000.000,  53.590.  245.200, 

*  9000.000,  -44.420,  9000.000, 

DATA  A10/  505.700,  9000.000, 

*  -340.200,  232.700,  9000.000, 

*  9000..  000,  9000.000,  304.100, 

*  9000.  000,  9000.000,  9000.000, 

*  751.900,  9000.000,  9000.000, 

*  9000.000,  9000.000,  9000.000, 

*  9000.000,  9000.000,  9000.000, 

*  9000. 000,  9000. 000,  9000. 000, 
DATA  All/  114.  SOO,  269.300, 

*  249.630,  10000.000,  -36.720, 

*  0.000,  372.900,  -235.700, 

*  9000. 000,  9000. 000,  9000. 000, 

*  9000.000,  108.900,  -209.700, 

*  9000.000,  9000.000,  183.000, 

*  -101.700,  148.300,  9000.000, 

*  9000.000,  -23.300,  9000.000, 

DATA  A12/  90.490,  91.650, 

*  155.700,  9000.000,  9000.000, 

*  -261.100,  0.000,  9000.000, 

*  9000.000,  9000.000,  9000.000, 

*  9000. 000,  9000. 000,  -287. 200, 

*  9000. 000,  9000. 000,  9000. 000, 

*  9000. 000,  9000. 000,  9000. 000, 

*  9000.000,  9000.000,  9000.000, 

DATA  A13/  83.360,  76.440, 

*  339.720,  -314.700,  9000.000, 

*  461.300,  9000.000,  0.000, 

*  9000.000,  9000.000,  9000.000, 

*  301.100,  137.800,  -154.300, 

*  95.180,  9000.000,  140.900, 

*  -20.110,  -149.500,  -202.300, 

*  9000.000,  145.600,  9000.000, 

DATA  A14/  -30.480,  79.400, 

*  -481.700,  -330.400,  9000.000, 

*  9000.000,  9000.000,  9000.000, 

*  -41.110,  9000.000,  9000.000, 

*  9000. 000,  9000. 000,  9000. 000, 

*  9000.000,  9000.000,  9000.000, 

*  9000.000,  9000.000,  9000.000, 

*  9000.000,  9000.000,  9000.000, 

DATA  A1 5/  65.330,  -41.320, 

*  -500.400,  -448.200,  9000.000, 

*  136.000,  9000.000,  -49.300, 

*  -189.200,  9000.000,  9000.000, 

*  9000. 000,  9000. 000,  9000. 000, 

*  9000.000,  9000.000,  9000.000, 

*  9000.000,  9000.000,  9000.000, 

*  9000.000,  9000.000,  9000.000, 

DATA  A16/  -83.980,  -186.000, 

*  -406.800,  -598.800,  9000.000, 


-46.280,  -163.700, 

-246. SOO,  -44. 580, 

-61.700,  9000.000/ 

9000. 000,  9000. 000,  -404. 800, 

128.000,  0.000, 

9000. 000,  9000. 000, 

9000. 000,  9000. 000, 

9000. 000,  9000. 000, 

9000. 000,  9000. 000, 

9000. 000,  9000, OCO, 

9000.000,  9000.000/ 

85.840,  -170.000,  245.400, 

372. 200,  9000. 000, 

9000. 000,  -73, 500, 

494. SOO,  660. 200, 

54.  570,  629.  000, 

9000. 000,  202. 300, 

9000. 000,  52. 080, 

9000. 000,  9000. 000/ 

9000.000,  9000.000,  191.200, 

9000. 000,  9000. 000, 

9000. 000,  SOOO. 000, 

9000. 000,  -356. 300, 

9000. 000,  9000. 000, 

4. 339,  9000. 000, 

9000. 000,  9000. 000, 

9000. 000,  9000. 000/ 

52.  130,  65. 690,  237. 700, 

52.  380,  -7.  838, 

9000. 000,  141. 700, 

9000. 000,  664. 600, 

47. 670,  9000. 000, 

-8. 538,  9000. 000, 

9000. 000,  172. 100, 

254.  SCO,  9000.000/ 

-44.850,  9000.000,  -164.000, 

9000. 000,  9000. 000, 

0.  000,  63.  720, 

9000. 000,  9000. 000, 

-99.810,  68.810, 

-70.140,  9000.000, 

9000. 000,  9000. 000, 

9000. 000,  9000. 000/ 

-22.310,  223.000,  -150.000, 

9000. 000,  9000. 000, 

ioa. aoo,  o. ooo, 

9000. 000,  9000. 000, 

71. 230,  4350. 000, 

9000. 000,  9000. 000, 

9000. 000,  9000. 000, 

9000. 000,  9000. 000/ 

-223.900,  109.900,  28.600, 

9000.  000,  9000. 000, 
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*  9000.000, 

*  0. 000, 

*  9000.000, 

*  9000. 000, 

*  9000.000, 

*  9000.000, 
DATA  A17/ 

*  5.  182, 

*  9000.000, 

*  9000.000, 

*  9000. 000, 

*  9000. 000, 

*  125.  300, 

*  9000. COO, 
DATA  A18/ 

*  -378. 200, 

*  9000. 000, 

*  9000. 000, 

*  9000. 000, 

*  9000. 000, 

*  9000. 000, 

*  9000. 000, 
DATA  A19/ 

*  157.800, 

*  -266. 600, 

*  9000.000, 

*  9000.000, 

*  9000.000, 

*  9000. 000, 

*  81.570, 
DATA  A20/ 

*  1020.000, 

*  -256. 300, 

*  9000. 000, 

*  44.420, 

*  9000.000, 

*  9000. 000, 

*  9000. 000, 
DATA  A21 / 

*  529.000, 

*  9000. 000, 

*  9000. 000, 

*  0.  000, 

*  9000.  000, 

*  9000. 000, 

*  9000. 000, 
DATA  A22/ 

*  669.  900, 

*  -132.900, 

*  -141.400, 

*  -84.530, 

*  9000.000, 

*  9000.000, 


9000.  000, 
9000. 000, 
-73. 850, 
9000. 000, 
9000. 000, 
9000. 000, 
5339.  000, 
-339. 500, 
9000. 000, 
0.  000, 
9000.  000, 
-62. 730, 
9000. 000, 
9000. 000, 
-101. 600, 
-332. 900, 
9000. 000, 
9000. 000, 
-351. 600, 
9000. 000, 
9000. 000, 
9000.  000, 
24.  820, 
242. 800, 
9000. OOO, 
617. 100, 
9000.000, 
9000. 000, 
9000. 000, 
-19. 140, 

315. 300, 
-66. 170, 
312. 500, 

9000. 000, 
-183. 400, 
9000. 000, 
9000. 000, 
-90.870, 
91. 460, 

698. 200, 
9000. 000, 
9000. 000, 

108. 300, 
9000. 000, 
9000. 000, 

-58. 770, 
34,  010, 
708. 700, 
9000. 000, 
9000. 000, 
0. 000, 
9000. 000, 
177. 600, 


9000.  000, 
9000.  000, 
-352.  900, 
9000.  000, 
9000. 000, 
9000. 000, 
9000. 000, 
9000. COO, 
9000. 000, 
9000. 000, 
9000. 000, 
9000. 000, 
9000. 000, 
9000. 000, 
9000.  000, 
-341.  600, 
9000.  000, 
0.  000, 
-114.  700, 
9000.  000, 
9000. 000, 
9000. 000, 
34.  780, 
9000. 000, 
9000.  000, 
134.  300, 
-15.  620, 
230.  900, 
9000. 000, 
9000.  000, 
349.  200, 
9000.  000, 
-338.  500, 
-313.  500, 
76.  750, 
9000.  000, 
9000. 000, 
9000.  000, 
-24.  360, 
9000.  000, 
225.  400, 
9000.  000, 
249.  200, 
450.  100, 
-125.  900, 
9000.  000, 
-52.  710, 
9000.  000, 
-197.  700, 
587.  300, 
0.  000, 
90G0. 000, 
9000,  QOO, 
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38.  890, 
9000. 000, 
-8.  283, 
9000.  000, 
9000. 000, 
9000. 000, 

650. 400, 
-399. 100, 
9000. 000, 
-216.  800, 
8455.  000, 
9000. 000, 
9000. 000, 
-293. 100, 
31. 870, 
-51. 540, 
9000. 000, 
-169. 700, 
-165. 100, 
9000. 000, 
9000. 000, 
9000. 000, 
-22. 970, 
-287. 500, 
9000. 000, 
0. 000, 
-54. 860, 
21. 370, 
-203. 000, 
9000. 000, 
62.  320, 
-297. 800, 
9000. OOQ, 
9000. 000, 
212. 700, 
9000. 000, 
9000. 000, 
9000. 000, 
4.  680, 

286. 300, 
9000. 000, 
9000. 000, 

62.  420, 
59. 020, 
9000. 000, 
9000. 000, 

121. 300, 

423. 200, 

9000. 000, 
9000. 000, 
56. 330, 
9000. 000, 
9000. 000, 


865. 900, 
9000. 000, 

-86. 360, 
9000. 000, 
9000. 000, 
9000. 000/ 

979. 800, 
9000. 000, 
9000. 000, 
9000. 000, 

699. 100, 
9000. 000, 
9000. 000, 
9000. 000/ 
49.  800, 
9000. 000, 
9000. 000, 
- 153.700, 
9000. 000, 
9000. 000, 
9000. 000, 
9000. 000/ 
-138. 400, 
9000. 000, 
9000. 000, 
9000. 000, 
52.  310, 
9000. 000, 
9000. 000, 
9000. 000/ 

268. 200, 
9000. 000, 
9000. 000, 
0. 000, 
9000. 000, 
9000. 000, 
9000. 000, 
9000. 000/ 

122. 900, 

-47, 510, 

9000. QOO, 

326. 400, 

464. 400, 
9000. 000, 
9000.  000, 
9000. 000/ 
9000. QOO, 
9000. 000, 
9000. 000, 
1821. 000, 
9000. 000, 
9000. 000, 

215. 000, 


529.  000, 


-132. 300, 


185.  400, 


-151. 000, 


562.  200, 


747.  700, 


♦ 

9000. 000, 

9000. 000, 

9000. 000, 

9000. 000, 

9000. 000/ 

DATA  A23/ 

36.  700, 

-185. 100, 

288. 500, 

33.  610, 

* 

649.  100, 

826. 760, 

9000. 000, 

552.  100, 

9000.  OOO, 

* 

176.  500, 

488. 900, 

-20. 930, 

9000. 000, 

9000.  000, 

* 

-293.  700, 

9000. 000, 

18.  980, 

74.  040, 

1346.  000, 

* 

-157.  100, 

0.  000, 

0.  000, 

-30.  100, 

9000.  000, 

* 

9000. 000, 

9000. 000, 

116. 600, 

9000. OOO, 

-64. 380, 

* 

9000. 000, 

86.  400, 

9000. 000, 

9000. 000, 

363. 700, 

* 

9000. 000, 

-79. 540, 

9000. 000, 

9000.  000, 

9000. 000/ 

DATA  A24/ 

-78. 450, 

-293.  700, 

-4. 700, 

134.  700, 

• 

860.  100, 

1201. 000, 

10000. 000, 

372.  000, 

9000.  000, 

* 

129.  500, 

9000. 000, 

113. 900, 

261.  100, 

91.  130, 

* 

-126.  000, 

1301. 000, 

309. 200, 

492.  000, 

689.  000, 

* 

ii.  aoo. 

17. 970, 

51. 900, 

0.  000, 

475.  800, 

• 

490.  900, 

534. 700, 

132. 200, 

9000. 000, 

546.  700, 

* 

9000. 000, 

247. 800, 

41.940, 

9000. 000, 

337.  700, 

• 

9000. 000, 

-86. 850, 

215. 200, 

498. 600, 

9000. 000/ 

DATA  A25/ 

-141. 300, 

-203. 200, 

-237. 700, 

375.  500, 

* 

661.  600, 

920. 400, 

9000. 000, 

128. 100, 

9000.  000, 

♦ 

-246.  300, 

9000. 000, 

9000. 000, 

203.  500, 

-108. 400, 

» 

ioaa.  ooo, 

323.  300, 

9000. 000, 

356. 900, 

9000. 000, 

• 

-314.  900, 

9000. 000, 

9000. 000, 

-255. 400, 

0.  000, 

«• 

-154.  500, 

9000. 000, 

9000. 000, 

9000. OOO, 

9000.  000, 

* 

9000.  000, 

9000. 000, 

-60. 700, 

9000. 000, 

9000. 000, 

9000.  000, 

9000. 000, 

9000.000, 

9000. 000, 

9000. 000/ 

DATA  A26/ 

-32.  690, 

-49.  920,. 

10.  380, 

-97.  050, 

♦ 

252.  600, 

417. 900, 

9000. 000, 

-142.  600, 

9000.  000, 

* 

9000.  000, 

9000. 000, 

-94. 490, 

9000. 000, 

9000.  000, 

* 

9000.  000, 

9000. 000, 

9000. 000, 

9000. OOO, 

9000. 000, 

* 

9000.  000, 

9000. 000, 

9000. 000, 

-34. 680, 

794.  400, 

♦ 

0.  000, 

533. 200, 

9000. 000, 

9000. OOO, 

9000. 000, 

* 

139.  800, 

304. 300, 

10. 170, 

-27. 700, 

9000.  000, 

• 

9000.  000, 

48.  400, 

9000. 000, 

9000. OOO, 

9000. 000/ 

DATA  A27 / 

5541. 000, 

9000. 000, 

1824. 000, 

-127.  800, 

W 

9000.  000, 

360. 700, 

9000. 000, 

9000. OOO, 

9000.  000, 

* 

9000.  000, 

9000. 000, 

9000. 000, 

9000. 000, 

9000.  000, 

* 

9000. 000, 

5250. 000, 

9000. 000, 

9000.  000, 

9000. 000, 

9000.  000, 

9000. 000, 

9000. 000, 

514. 600, 

9000.  000, 

♦ 

-85. 120, 

0.  000, 

9000. 000, 

9000. 000, 

9000.  000, 

♦ 

9000.  000, 

9000. 000, 

9000. 000, 

9000.  000, 

9000. 000, 

* 

9000. 000, 

9000. 000, 

9000, 000, 

9000. OOO, 

9000. 000/ 

DATA  A28/ 

-52. 650, 

16. 620, 

21. 500, 

40.  680, 

» 

914.  200, 

1081. 000, 

9000. 000, 

303. 700, 

9000. 000, 

* 

243.  800, 

9000.  000, 

112.  400, 

9000. 000, 

9000.  000, 

* 

9000. 000, 

9000. 000, 

9000. 000, 

335.  700, 

9000.  000, 

* 

-73.  090, 

9000. 000, 

-26. 060, 

-60. 710, 

9000.  000, 

♦ 

9000.  000, 

9000. 000, 

0.  000, 

9C00.  000, 

9000.  000, 

* 

9000. 000, 

9000. 000, 

9000. 000, 

9000. OOO, 

9000.  000, 

» 

9000. GGG, 

a  •**»  r\ 

~  *±  /  .  J  t 

nnno  npn 
JUUU • uuu, 

3UUU • UUUf 

npnp  f~s / 

2UUU.  UllU/ 

DATA  A29/ 

-7.  481, 

9000. 000, 

28. 410, 

9000. 000, 

# 

382. 800, 

9000.  000, 

9000. 000, 

1 60 . 600, 

9000. 000, 

# 

9000. 000, 

239.  800, 

63.  710, 

106. 700, 

9000. 000, 

♦ 

9000.  000, 

9000.  000, 

9000. 000, 

125.  700, 

9000. 000, 

742.  ICO, 


856. 300, 


246.  900, 


341. 700, 


561. 600, 


823. 500, 


461. 600, 
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*  -27. 940, 

*  9000.000, 

*  9000.000, 

*  9000.000, 
DATA  A3G/ 

*  9000.000, 

*  -146.300, 

*  9000.000, 

*  9000.000, 

*  9000. 000, 

*  9000.000, 

*  9000. 000, 
DATA  A31/ 

*  9000.000, 

*  152. 000, 

*  9000.000, 

*  9000. 000, 

*  481.300, 

*  0.000, 

*  9000.000, 
DATA  A32/ 

*  9000. 000, 

*  21.920, 

*  9000. 000, 

*  9000.000, 

*  64. 280, 

*  9000.000, 

*  9000.000, 
DATA  A33/ 

*  9000. 000, 

*  9000.000, 

*  9000.000, 

*  1169.000, 

*  125.300, 

*  9000. 000, 

*  9000.000, 
DATA  A34/ 

*  9000.000, 

*  9000. 000, 

*  9000. 000, 

*  9000.000, 

*  174.400, 

*  9000. 000, 

*  9000. 000, 
DATA  A35/ 

*  695.000, 

*  41.570, 

*  9000.000, 

*  9000.000, 

*  9000.000, 

*  535.800, 

*  9000. 000, 
DATA  A36/ 


9000. 000, 
9000. 000, 
9000. 000, 
9000. 000, 
-25.  310, 
23.  480, 
9000. 000, 
9000. 000, 
9000. 000, 
9000. 000, 
9000. 000, 
9000. 000, 
140. 000, 
0. 000, 
9000. 000, 

164. 400, 
9000. 000, 
9000. 000, 
9000. 000, 
9000. 000, 
128. 000, 
9000. 000, 
9000. 000, 
9000. 000, 
-40. 820, 
9000. 000, 
0.  000, 
9000. 000, 
-31. 520, 
9000. 000, 
9000. 000, 
9000. 000, 
9000.  000, 
9000. 000, 
9000. 000, 
9000. 000, 
-72.  880, 
9000. 000, 
9000.  000, 
9000. 000, 
9000. 000, 
9000. 000, 
9000.  000, 
9000. 000, 
50. 490, 
-240. 000, 
9000. 000, 
9000. 000, 
-215.  000, 
9000. 000, 
9000. 000, 
9000. 000, 
-165. 900, 


9000. 000, 
9000. 000, 
9000. 000, 
9000. 000, 
9000. 000, 
9000. 000, 
9000. 000, 
9000. 000, 
48. 480, 
9000. 000, 
9000. 000, 
9000. 000, 
9000. 000, 

838. 400, 
9.  207, 
9000. 000, 
9000. 000, 
9000. 000, 
9000. 000, 
9000. 000, 
9000.000, 
9000. 000, 
476.  600, 
9000. 000, 
21. 760, 
9000.000, 
9000. 000, 
9000.000, 
9000. 000, 
9000. 000, 
736.  400, 
9000.  000, 
9000. 000, 
9000. 000, 
0.  000, 
9000.  000, 
-184. 400, 
9000. 000, 
9000. 000, 
9000. 000, 
9000. 000, 
9000, 000, 
9000. 000, 
9000. 000, 
9000. 000, 
9000. 000, 
-122. 100, 
9000.  000, 
-343. 600, 
9000. 000, 
9000. 000, 
9000. 000, 
9000.  000, 


9000. 000, 
0.  000, 
9000. 000, 
78. 920, 
157. 300, 
317. 500, 
9000. 000, 
9000.  000, 
-133.  160, 
9000.  000, 
9000. 000, 
9000. 000, 
221.  400, 
9000. 000, 
9000.  OOO, 
9000. 000, 
9000. 000, 
9000. 000, 
9000. 000, 
302. 200, 
58.  680, 
138. 000, 
9000. 000, 
9000. 000, 
48. 490, 
9000. 000, 
9000. OOO, 
9000. OOO, 
155.  600, 
-142. 600, 
9000. OOO, 
9000. OOO, 
225.  800, 
9000. OOO, 
9000. OOO, 
9000. OOO, 
9000. OOO, 
443.  600, 
9000. OOO, 
329. 100, 
9000. OOO, 
9000. OOO, 
0. OOO, 
-119. 300, 
-2.  504, 

110. 400, 
9000. OOO, 
9000. OOO, 
-58. 430, 
85. 700, 
9000. OOO, 
-97. 710, 
9000. OOO, 


9000. OOO, 

9000. 000, 

31. 660, 

9000. 000/ 

404. 300,  521 

9000. OOO, 

9000. OOO, 

9000. OOO, 

9000. OOO, 

0. OOO, 

9000. OOO, 

9000. 000/ 

150. 600,  267 

9000. OOO, 

9000. OOO, 

9000. OOO, 

9000. OOO, 

9000. OOO, 

-417. 200, 

9000. 000/ 
9000.000,  501 

9000. OOO, 

9000. OOO, 

9000. OOO, 

9000. OOO, 

9000. OOO, 

9000. OOO, 

9000. 000/ 

291.  lOO,  721. 
9000. OOO, 

9000. OOO, 

9000. OOO, 

224. OOO, 

9000. OOO, 

9000. OOO, 

9000. 000/ 

9000. OOO,  9000. 
9000. OOO, 

9000. OOO, 

9000. OOO, 

9000. OOO, 

9000. OOO, 

9000. OOO, 

9000. 000/ 
-143.200,  -25. 

9000. OOO, 

9000. OOO, 

9000. OOO, 

9000. OOO, 

9000. OOO, 

0. OOO, 

9000. 000/ 

9000. OOO,  9000. 


.  600, 


600, 


300, 


900, 


OOO, 


870, 


OOO, 
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*  9000.000,  386.600,  9000.000, 

*  9000.000,  9000.000,  9000.000, 

*  9000.000,  9000.000,  9000.000, 

*  9000.000,  9000.000,  9000.000, 

*  9000.000,  9000.000,  9000.000, 

*  9000. 000,  9000. 000,  9000. 000, 

*  0. 000,  9000. 000,  9000. 000, 

DATA  A37/  41.900,  -3.167, 

*  726.700,  9000.000,  9000.000, 

*  -34.140,  9000.000,  -209.300, 

*  9000.000,  9000.000,  9000.000, 

*  201.700,  9000.000,  85.320, 

*  313.800,  9000.000,  167.900, 

*  9000.000,  9000.000,  9000.000, 

*  9000.000,  0.000,  9000.000, 

DATA  A38/  -5.132,  9000.000, 

*  645.900,  9000.000,  9000.000, 

*  9000.000,  9000.000,  9000.000, 

*  9000.000,  9000.000,  9000.000, 

*  9000. 000,  9000. 000,  9000. 000, 

*  9000. 000,  9000. 000,  9000. 000, 

*  9000.000,  9000.000,  9000.000, 

*  9000.  000,  9000. 000,  0. 000, 

DATA  A39/  -31.950,  37.700, 

*  172.  2C0,  -287.100,  9000.000, 

*  9000.000,  9000.000,  -158.200, 

*  9000.000,  335.600,  9000.000, 

*  9000. 000,  9000. 000,  9000. 000, 

*  9000. 000,  9000. 000,  900p. 000, 

*  -191.700,  9000.000,  9000.000, 

*  9000.000,  9000.000,  9000.000, 

DATA  A40/  147.300,  9000.000, 

*  9000.000,  9000.000,  9000.000, 

*  9000.  COO,  9000.  000,  9000.  000, 

*  9000.000,  9000.000,  9000.000, 

*  9000.000,  9000.000,  9000.000, 

*  9000.000,  9000.000,  9000.000, 

*  9000.000,  9000.000,  9000.000, 

*  9000.000,  9000.000,  9000.000, 
DO  5  1=1, 40 

AKI,  1)=A1(I) 

AI ( I,  2) =A2 ( I  ) 

AI ( I,  3) =A3< I ) 

AKI,  4)  =  A4  ( I  ) 

Aid,  5)  =A5(I) 

AI ( I,  6) =A6<I  ) 

AI (I,  7 ) = A7 ( I  ) 

AI d,  8) =A8< I  ) 

AI  (I,  9  )  =  A9  <  I  ) 

AI  (  I,  10  )  =  A 1 0  <  I  ) 

Aid,  11  )  =  A 1 1  (  I  ) 

AI C I,  12) =A12 ( I  ) 

AI C I,  13) *A13< I  ) 


9000.000,  9000.000, 

9000.000,  9000.000, 

-42.310,  9000.000, 

9000.000,  9000.000, 

9000.000,  9000.000, 

9000.000,  9000.000, 

9000.000,  9000.000/ 

-75.670,  9000.000,  640.900, 

-8.671,  9000.000, 

9000.000,  9000.000, 

298.400,  2344.000, 

143.200,  9000.000, 

9000.000,  9000.000, 

9000.000,  9000.000, 

9000.000,  9000.000/ 

-237.200,  -157.300,  649.700, 

9000.000,  9000.000, 

9000.000,  9000.000, 

9000.000,  9000.000, 

-124.600,  9000.000, 

9000.000,  9000.000, 

9000.000,  9000.000, 

9000.000,  9000.000/ 

-133.900,  -240.200,  64.160, 

97.040,  9000.000, 

9000. 000,  9000. 000, 

9000. 000,  9000. 000, 

-186.700,  9000.000, 

-71.000,  9000.000, 

6.699,  136.600, 

0.000,  9000.000/ 

9000.000,  9000.000,  9000.000, 
9000.000,  9000.000, 

9000.000,  9000.000, 

9000.000,  9000.000, 

9000.000,  9000.000, 

9000.000,  9000.000, 

9000.000,  9000.000, 

9000.000,  0.000/ 
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c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


AKI,  14  )  “A14  <  I ) 
AKI,  15)  =A15(I ) 

A I  <  I,  16)  =A16d  ) 
AI<I, 17)=A17(I) 
Aid,  ia)=Aia<i> 
AKI,19)=A19d> 
AKI,  20)=A20d) 
AKI,  21)*A21d) 
AKI,  22)=A22d) 
Aid,  23 )  =A23  ( I  > 
AKI,  24  )  -A24  ( I ) 
AKI,  25)  =A25d) 
Aid,  26)=A26d> 
■AKI,27)«A27d> 

AI < I, 2S ) =  A28 ( I ) 
Aid,  29)  =A29d> 
Aid,  30 ) =A30 ( I ) 
AKI,  31)  =A31d> 
AKI,  32 )  aA32  ( I ) 
AI<I, 33 ) ® A33  < I ) 
Aid,  34)  =A34<I) 

AI ( I,  35) =A35(I ) 
AKI,  36 )  =A36  <  I ) 

AI < I, 37 ) =A37 ( I ) 
AKI,3a>=A38d) 
AKI,  39>=A39<I) 

AI ( I,  40 ) aA40 ( I ) 

5  CONTINUE 

DO  10  1  =  1,76 
RK  I  )  =RR  <  I  ) 

QI ( I ) =QQ ( I > 

MAlNd  )  =«AINSG(I ) 
10  CONTINUE 
RETURN 
END 


ULLE  CONTAINS  BUILT-IN  UNIFAC  LLE- PARAMETERS 


* 


SUBROUTINE  ULLE < RI, QI, AI, MAIN > 


C  . . . . 

c 

C  THE  MAIN  GROUPS  ARE: 


c 

1 

CH2  . . . 

.  .  .  .  2 

C  =  C . 

.  3 

ACH . . 

.  4 

ACCH2  . . . 

.  5 

OH  .  .  . 

c 

6 

PI  .... 

7 

P2  ...... 

.  a 

H20 . . 

.  9 

ACOH  . . . . 

10 

CH2C0 

c 

11 

CHO  . . . 

.  .  .  12 

FURFURAL 

13 

COOH  .... 

14 

COOC  . . . . 

15 

CH2Q  . 
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C  16  CCL .  17 

C  21  CCN .  22 

C  26  DEOH .  27 

C  31  TMS .  32 

r 


CCL2  .... 

18 

CCL3 

ACNH2  . . . 

23 

CN02 

PYRIDINE 

28 

TCE  . 

DMSO  .... 

19 

CCL4  .... 

20 

ACCL 

24 

ACNQ2  . . . 

25 

DOH  . 

29 

MFA  . 

30 

DMF  A 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


THE  SUB  GROUPS  ARE: 


1 

CH3 . 

2 

CH2 . 

2 

CH . 

4 

C . 

.  5 

CH2=CH  . . . 

8 

CH=CH  .... 

7 

CH=C  .... 

a 

CH2=C  . . . . 

9 

ach . 

10 

AC . 

11 

ACCH3  .... 

12 

ACCH2  . . . 

13 

ACCH  .... 

14 

OH . 

15 

PI . 

16 

P2 . 

17 

H20 . 

18 

ACOH  .... 

19 

CH3C0  . . . 

20 

CH2C0  . . . . 

21 

CHO . 

22 

FURFURAL 

23 

COOH  .... 

24 

HCOOH  . . . 

25 

CH3C00  . . . 

26 

CH2C00  . . . 

27 

CH30  .... 

28 

CH20  .... 

29 

CH-0 - 

30 

FCH20  . . . . 

31 

CH2CL  .... 

32 

CHCL _ 

33 

CCL  « •  •  •  • 

34 

CH2CL2  . . 

35 

CHCL2  _ 

36 

CCL.2 . 

37 

CHCL3  . . . 

38 

CCL3  .... 

39 

CCL4  .... 

40 

ACCL . 

41 

CH3CN  .... 

42 

CH2CN  . . . 

43 

ACNH2  . . . 

44 

CH3N02  . . 

45 

CH2N02  . . . 

46 

CHN02  .... 

47 

ACN02  . . . 

48 

( CH20H ) 2 

49 

<  HOMM ) 20 

50 

C5H5N  . . . . 

51 

C5H4N  .... 

52 

C5H3N  . . . 

53 

CCL2=CHCL 

54 

HCQNHCH3 

55 

HCON ( CH3 ) 1 

56 

( CH2 ) 4S02 

57 

( CH2 ) 2S0 

IMPLICIT  REAL*8< A-H. Q-Z) 

DIMENSION  AX  <  40, 40  > , RR  <  57 ) , QQ  <  57 ) ,  RI ( 76 ) ,  QI < 76 > , MAINSG ( 57 ) , 

*  MAIN (76) 

DIMENSION  Al(32),  A2<32),  A3 (32),  A4<32>,  A5<32),  AS (32),  A7 ( 32  > , Afi  <  32 ) 

*  A9<  32 ) , A10  <  32 ) , All <  32 ) ,  A12  <  32 ) ,  A13  ( 32 ) ,  A14  ( 32 > ,  A15 < 32 ) , A16 ( 32 ) , 
*A17 ( 32 ) ,  A18<32), AI9<32),  A20C32),  A21 < 32 ) , A22 < 32 > ,  A23 ( 32 ) ,  A24 ( 32 ) , 
*A25<32), A26 ( 32 ) , A27 ( 32 ) ,  A28 ( 32 ) ,  A29 ( 32 ) , A30 ( 32 > ,  A31(32>,  A32(32) 

DATA  MAINSG/  4*1,  4*2,  2*3,  3*4,  1*5,  1*6,  1*7,  1*8,  1*9,2*10, 

*  1*11,  1*12,  2*13,  2*14,  4*15,  3*16,  3*17,  2*  18,  1*19,  1*20, 

*  2*21,  1*22,  3*23,  1*24,  1*25,  1*26,  3*27,  1*2Q,  1*29,  1*30, 


» 


*  1*31, 1*32/ 

.DATA  RR/G.  9011,  0.  6744,  0,  4469,  0.  2195,  1.  3454,  1.  1167,  0.  8886,  1.  1173, 
*0.  5313,  0.  3652,  1.  2663,  1.  0396,  0.  8121,  1. 0000,  3.  2499,  3,  2491,  0.  9200, 
*0.  8952,  1.  6724,  1.  4457,  0.  9980,  3.  1680,  1.  3013,  1.  5280,  1.  9031,  1. 6764, 

*1,  1450,  O.  9183,  0.  6908,  0.  9183,  1.  4654,  1.  2380,  1.  0060,  2.  2564,  2.  0606, 

*1. 8016,  2.  8700,  2.  6401,  3.  3900,  1.  1562,  1.  8701,  1.  6434,  1.  0600,  2.  0086, 

*1.  7Q18,  1.  5544,  1.  4199,  2.  4088,  4.  0013,  2.  9993,  2.  8332,  2.  6670,  3.  3092, 

*2.  4317,  3.  0856,  4.  0358,  2.  8266/ 

DATA  QQ/O.  848,  0.  540,  0.  228,  0.  000,  1.  176,  O.  867,  0.  676,  0.  988,  0.  400, 
*0.  120,  0.  968,  0.  660,  0.  348,  1. 200,  3.  128,  3.  124,  1.  400,  0.  680,  1.  488, 

*1.  180,  0.  948,  2.  484,  1.  224,  1. 532,  1. 728,  1.  420,  1.  088,  0.  780,  0.  468, 

*1.  100,  1.  264,  0.  952,  0.  724,  1.  988,  1. 684,  1.  448,  2.  410,  2.  184,  2.  910, 

*0.  844,  1.  724,  1.  416,  0.  816,  1. 868,  1.  560,  1.  248,  1.  104,  2,  248,  3,  568, 

*2.  113,  1.  833,  1.  553,  2.  860,  2.  192,  2.  736,  3.  200,  2.  472/ 

DATA  Al/  0.000,  74.540,  -114.800,  -115.700,  644.600, 

*  329.600,  310.700,  1300.000,  2255.000,  472.600, 

972. 400, 


158. 100, 
42.  140, 
696. 800, 
269. 200, 
561. 400, 


383.  000, 
■243. 900, 
902.  200, 
•300. 000, 
956. 500/ 


139. 400, 
7. 500, 
556. 700, 
-63. 600, 


-5.  550, 
575.  700, 
928.  300, 


662.  100, 
924.  800, 
527.  500, 
331. 000, 
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DATA  A2/ 

292. 300, 

0. 000, 

340.  700, 

4102. 000, 

724. 400 

1731. 000, 

1731. 000, 

396. 000, 

9000.  000, 

343. 700, 

-214. 700, 

9000. 000, 

1647.  000, 

-577. 500, 

289. 300, 

99.  610, 

337. 100, 

4583. 000, 

5831.  000, 

9000. 000, 

405.  900, 

9000. 000, 

425.  700, 

9000.  000, 

9000. 000, 

9000’.  000, 

9000. 000, 

9000.  000, 

500.  700, 

115.  400, 

784.  400, 

265. 400/ 

DATA  A3/ 

156. 500, 

-94. 780, 

0.  000, 

167. 000, 

703. 900 

511.  500, 

577. 300, 

859.  400, 

1649.  000, 

593. 700, 

362. 300, 

31. 140, 

461.  800, 

6.  000, 

32.  140, 

-is. aio. 

9000.  000, 

-231. 900, 

3000. 000, 

-878. 100, 

29.  130, 

1.  640, 

-1. 770, 

-11. 190, 

358. 900, 

363.  500, 

-573.  200, 

9000.  000, 

364.  200, 

-58.  100, 

21.  970, 

84. 160/ 

DATA  A4/ 

104.  400, 

-269.  700, 

-146. 600, 

0.  000, 

4000.  000 

136. 600, 

906. 800, 

5695.  000, 

292.  600, 

916. 700, 

1218. 000, 

715, 600, 

339. 100, 

5688.  000, 

213. 100, 

-114. 100, 

9000. 000, 

-12. 140, 

-141.  300, 

-107. 300, 

1208. 000, 

689. 600, 

3629.  000, 

-175.  600, 

337. 700, 

1023. 000, 

-390. 700, 

9000. 000, 

9000. 000, 

9000. 000, 

238. 000, 

132. 200/ 

DATA  A5/ 

328. 200, 

470.  700, 

-9.  210, 

1.  270, 

0.  000 

937. 300, 

991.  300, 

28.  730, 

-195. 500, 

67. 070, 

1409. OCO, 

-140. 300, 

-104.  000, 

195.  600, 

262.  500, 

62.  050, 

272. 200, 

-61.570, 

-41.  750, 

-597. 100, 

-189. 300, 

-348. 200, 

-30. 700, 

-159.  000, 

536. 600, 

53.  370, 

183.  300, 

-44.  440, 

9000. 000, 

9000. 000, 

9000. 000, 

9000. 000/ 

DATA  A6/ 

-136. 700, 

-135.  700, 

-223.  000, 

-162. 600, 

-281. 100 

0. 000, 

0.  000, 

-61. 290, 

-153.  200, 

-47. aio. 

-344. 100, 

299.  300, 

244. 400, 

19.  570, 

1970. 000, 

-166. 400, 

128.  600, 

1544. 000, 

224.  600, 

9000. 000, 

9000. 000, 

9000. 000, 

150.  800, 

9000. OGO, 

9000. 000, 

9000. 000, 

9000. OCO, 

9000. 000, 

9000. 000, 

9000. 000, 

9000. 000, 

9000. 000/ 

DATA  A7 / 

-131.  900, 

-135.  700, 

-252. 000, 

-273. 600, 

-268. 800 

0. 000, 

0.  000, 

5.  890, 

-153. 200, 

353. 800, 

-338. 600, 

-241.  800, 

-57. 980, 

487.  100, 

1970. 000, 

-166. 400, 

507.  800, 

1544. 000, 

-207. 000, 

9000. 000, 

9000. 000, 

9000.  000, 

150. 800, 

9000. 000, 

9000. 000, 

9000. 000, 

9000. 000, 

9000. 000, 

9000. 000, 

9000. 000, 

9000. 000, 

9000.  000/ 

DATA  AS/ 

342.  400, 

220. 600, 

372. 800, 

203. 700, 

-122. 400 

247. 000, 

104. 900, 

0.  000, 

344.  500, 

-171. 800, 

-349. 900, 

66.  950, 

-465.  700, 

-6.  320, 

64.  420, 

315. 900, 

370. 700, 

356. 800, 

502. 900, 

-97.  270, 

193.  300, 

-109. 800, 

1539. 000, 

32.  920, 

-269. 200, 

9000. 000, 

-87 J.  600, 

1429.  000, 

-364.  200, 

-117. 400, 

13. 410, 

9000.  000/ 

DATA  A9/ 

-159.  800, 

9000. 000, 

-473. 200, 

-470. 400, 

-63.  150 

-547. 200, 

-547.  200, 

-595. 900, 

0.  000, 

-825.  700, 

9000. 000, 

9000.  000, 

9000.  000, 

-898. 300, 

9000. 000, 

9000. 000, 

9000. 000, 

9000.  000, 

4894.  000, 

9000. 000, 
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*  9000.000,  -851.600, 

*  9000.000,  -S37. 300. 

*  9000.000,  9000.000/ 

DATA  A  IQ/  66.560, 

*  401.700,  -127.600, 

*  -37.360,  120.300, 

*  1000.000,  -301.000, 

*  430.600,  1010.000, 

*  9000.000,  9000.000, 

*  9000.000,  9000.000/ 

DATA  All/  146. 100, 

*  643.400,  231.400, 

*  0.000,  1724.000, 

*  751.800,  9000.000, 

*  9000. 000,  9000. 000, 

*  9000. 000,  9000. 000, 

*  9000.000,  9000.000/ 

DATA  A12/  14. 780, 

*  -94.640,  732.300, 

*  -311.600,  0.000, 

*  9000. 000,  -347. 900, 

*  9000. 000,  9000. 000, 

*  9000. 000,  -208. 400, 

*  9000. 000,  9000. 000/ 

DATA  A13/  1744.000, 

*  728.  700,  349.  100, 

*  1051.000,  -115.700, 

*  19.770,  1670.000, 

*  9000. 000,  942. 200, 

*  9000. 000,  9000. 000, 

*  9000.000,  9000.000/ 

DATA  A14/  -320. 100, 

*  -76.640,  -152.800, 

*  1090.000,  -46.130, 

*  9000.000,  108.900, 

*  -149.200,  9000.000, 

*  9000. 000,  9000. 000, 

*  9000. 000,  9000. 000/ 

DATA  A15/  1571.000, 

*  -218. 100,  -218. 100, 

*  9000. 000,  9000. 000, 

*  301.100,  137.800, 

*  9000.000,  9000.000, 

*  9000. 000,  9000. 000, 

*  9000. 000,  9000. 000/ 

DATA  A16/  73. 800, 

*  351.500,  351.500, 

*  -47.510,  9000.000, 

*  0.000,  110.500, 

*  9000.000,  9000.000, 

«  9000.000,  9000.000, 

*  9000.000,  9000.000/ 

DATA  A 17/  27.900, 


9000.000,  -16.130, 

9000. 000,  9000. 000, 

306. 100,  -78. 310, 

634.800,  -568.000, 

1247.000,  258.700, 

12.010,  -10.880, 
400.000,  -328.600, 

148. 000,  9000. 000, 

517.000,  -75.300, 

623. 700,  9000. 000, 

0.750,  -245.800, 

9000.000,  9000.000, 
9000. 000,  9000. 000, 
9000.000,  9000.000, 

9000.000,  -10.440, 

211.600,  9000.000, 

1919.000,  57.700, 

-249.300,  61.590, 

9000. 000,  9000. 000, 

-13.910,  90Q0.  000, 

-48.520,  75.490, 

652. 300,  9000. OOO, 
0.000,  -117.600, 

48. 150,  43.  830, 

446. 300,  9000. OOO, 
-2.160,  9000.000, 

485.600,  114.800, 

385.900,  -337.300, 

1417.000,  0.000, 

-209.700,  54.570, 

9000.000,  9000.000, 
9000. OOO,  9000. OOO, 

76. 440,  52. 130, 

212.800,  9000.000, 

1402.000,  461.300, 

-154.300,  47.670, 

95.180,  9000.000, 
9000.000,  9000.  OOO, 

-24.  360,  4.  680, 

770.000,  9000.000, 

337.100,  9000.000, 

249.200,  62.420, 

9000.000,  9000.000, 

9000.000,  9000.000, 

-52.710,  9000.000, 


-538.  600, 

9000. OOO, 

-73.870,  216.000, 

0.  OOO, 

5.  202, 

902.  600, 

211.  600, 

9000.  OOO, 

223.200,  -431.300, 

12S.  OOO, 

9000. OOO, 

9000.  OOO, 

9000.  OOO, 

9000.  OOO, 

-184.900,  444.700, 

48.  930, 

9000. OOO, 

9000. OOO, 

-278.  200, 

173.  800, 

147.300,  118.400, 

-101.  300, 

-96.  620, 

874.  300, 

572.  700, 

9000.  OOO, 

-170.000,  180.600, 

58. 840, 

-235.  700, 

629.  OOO, 

343.  100, 

9000. OOO, 

65.690,  137.100, 

52.  380, 

0.  OOO, 

9000. OOO, 

9000. OOO, 

9000. OOO, 

122.900,  455.100, 

483.  900, 

225.  400, 

9000.  OOO, 

9000.  OOO, 

9000.  OOO, 

9000. OOO,  669. 200, 
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♦ 

-ias.  loo, 

-401. 600, 

740. 400, 

9000. 000, 

550. 600, 

♦ 

9000.  000, 

aos. aoo. 

437. 700, 

-132. 900, 

-197. 700, 

• 

-21.  350, 

0.  000, 

9000. QOO, 

56. 330, 

9000. 000, 

• 

9000. 000, 

9000. 000, 

9000. 000, 

9000. 000, 

9000. 000, 

* 

9000. 000, 

9000.  000, 

9000. 000, 

9000. 000, 

9000. 000, 

♦ 

9000'.  000, 

9000. 000/ 

DATA  A18/ 

21.  230, 

-185. 100, 

288. 500, 

33.  610, 

418.  400, 

# 

-465.  700, 

-465. 700, 

793. 200, 

9000. 000, 

342. 200, 

* 

9000. 000, 

203. 100, 

370. 400, 

176. 500, 

-20.  930, 

* 

-157.  100, 

9000. 000, 

0. 000, 

-30. 100, 

9000.  000, 

» 

70.  040, 

-75.  500, 

9000. 000, 

9000. 000, 

9000. 000, 

* 

9000. 000, 

18.  980, 

9000. 000, 

9QQ0. 000, 

9000. 000, 

*• 

9000. 000, 

9000. 000/ 

DATA  A19/ 

39.  970, 

-293. 700, 

-4.  700, 

134.  700, 

713.  500, 

* 

-260. 300, 

512.  200, 

1205.  000, 

1616. 000, 

550.  000, 

♦ 

9000.  000, 

70.  140, 

433. 100, 

129. 500, 

113.  900, 

♦ 

li.  aoo. 

17. 970, 

51. 900, 

0.  000, 

475.  800, 

♦ 

492.  000, 

1302. 000, 

490. 900, 

534. 700, 

9000. 000, 

* 

9000. 000, 

9000. 000, 

9000. 000, 

9000. 000, 

9000. 000, 

• 

9000. 000, 

9000. 000/ 

DATA  A20/ 

-59.  060, 

9000. 000, 

777. 800, 

-47.  130, 

1989. OOO, 

* 

9000. 000, 

9000. 000, 

390. 700, 

9000. 000, 

190.  500, 

» 

9000. 000, 

9000. 000, 

1349. 000, 

-246. 300, 

9000.  000, 

♦ 

9000. 000, 

9000. QOO, 

9000. 000, 

-255. 400, 

0.  000, 

# 

346.  200, 

9000. 000, 

-154. 500, 

9000. 000, 

124.  aoo. 

* 

9000.  000, 

-387. 700, 

9000. 000, 

9000, 000, 

9000. OOO, 

* 

9000. 000, 

9000. 000/ 

DATA  A21 / 

29. 080, 

34. 780, 

56. 410, 

-53.  290, 

2011.  OOO, 

* 

9000.  000, 

9000. QOO, 

S3. 480, 

9000. OCO, 

-349.  200, 

i 

• 

9000. 000, 

9000. 000, 

9000. 000, 

2.  410, 

9000. OOO, 

• 

9000. OOO, 

9000. 000, 

-15. 620, 

-54. 860, 

-465. 200, 

♦ 

0.  000, 

9000. 000, 

9000. 000, 

9000. QOO, 

9000. OOO, 

♦ 

9000. 000, 

134. 300, 

9000. 000, 

9000. QOO, 

9000. OOO, 

# 

9000.  000, 

9000. 000/ 

DATA  A22/ 

175.  800, 

9000. 000, 

-218. 900, 

-15.  410, 

529.  OOO, 

# 

9000. 000, 

9000.  000, 

-239. aoo. 

-860. 300, 

857.  700, 

• 

9000.  000, 

9000. 000, 

681. 400, 

9000. 000, 

9000. OOO, 

• 

9000. 000, 

9000. 000, 

-216. 300, 

8455. 000, 

9000.  OOO, 

9000. 000, 

0.  000, 

9000. 000, 

179. 900, 

125.  300, 

* 

9000. 000, 

924.  500, 

9000. 000, 

9000. 000, 

9000. OOO, 

9000. 000, 

9000. 000/ 

DATA  A23/ 

94.  340, 

375. 400, 

113. 600, 

-97. 050, 

483.  800, 

■* 

264.  700, 

264. 700, 

13.  320, 

9000. 000, 

377. OOO, 

♦ 

9000. 000, 

9000. 000, 

152.  400, 

9000. 000, 

-94.  490, 

* 

9000.  000, 

9000. 000. 

9000. 000, 

-34. 680, 

794. 400, 

* 

9000.  000, 

9000. 000, 

0. 000, 

9000. 000, 

139.  800, 

* 

9000. 000, 

9000. 000, 

9000. 000, 

9000. 000, 

9000. OOO, 

* 

9000. 000, 

9000. 000/ 

DATA  A24/ 

193.  600, 

9000. 000, 

7. 180, 

-127. 100, 

332.  600, 

* 

9000. 000, 

9000. 000, 

439. 900, 

-230. 400, 

211. 600, 

* 

9000. 000, 

9000. 000, 

9000. 000, 

9000. 000, 

9000. OOO, 

♦ 

9000.  000, 

9000. 000, 

9000. 000, 

514. 600, 

9000. OOO, 

• 

9000.  000, 

175.  aoo. 

9000. 000, 

0.  000, 

963. OOO, 
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9000. 000, 
9000. 000, 
DATA  A25 / 
9000. 000, 
9000. 000, 
9000. 000, 
9000. 000, 
9000. 000, 
9000. 000, 
DATA  A26/ 
9000. 000, 
9000. 000, 
9000. 000, 
9000.  000, 
0.  000, 
9000.  000, 
DATA  A27 / 
9000. 000, 
9000.  000, 
9000. 000, 
-169. 700, 
9000. 000, 
9000. 000, 
DATA  A2S/ 
9000. 000, 
9000. OOO, 
9000. 000, 
9000. 000, 
£000.  000, 
9000. 000, 
DATA  A29/ 
9000. 000, 
9000. 000, 
9000. 000, 
9000. 000, 
9000. 000, 
9000, QOQ, 
DATA  A30/ 
9000. 000, 
9000. 000, 
9000. 000, 
9000. 000, 
9000. 000, 
9000. 000, 
DATA  A31 / 
9000. 000, 
9000. 000. 
9000. 000, 
9000. 000, 
9000. 000, 
0. 000, 
DATA  A32 / 
9000.  000, 


9000. 000, 
9COO. 000/ 
108.  500, 
9000. 000, 
-75. 230, 
9000.  000, 
164.  400, 
9000.  000, 
9000. 000/ 
Si.  490, 
9000. 000, 
9000. 000, 
9000.  000, 
9000. 000, 
9000.000, 
9000. 000/ 
-128. 800, 
9000. 000, 
-201. 900, 
9000. 000, 
-944. 900, 
0.  000, 
9000. 000/ 
147.  300, 
9000.000, 
123. 200, 
9000. 000, 
9000.000, 
9000.  000, 
9000. OOO/ 
-11. 910, 
9000. OOO, 
9000. OOO, 
9000. OOO, 
9000.  OOO, 
9000.  OOC, 
9000. 000/ 
14.  910, 
9000. OOO, 
-281.  900, 
9000. OOO, 
9000. 000, 
9000.  OOO, 
9000. 000/ 
£7.  S40, 
9000.  OOO, 
9000.  OOO, 
9000.  000, 
9000. OOO, 
9000. OOO, 
9000. 000/ 
36.  420, 
9000. OOO, 


9000. OOO, 

9000. OOO, 
-424.  300, 
1707. OOO, 
9000. 000, 
481.  300, 
9000.  000, 

90Q0.  OOO, 
9000.  OOO, 
9000.  OOO, 
9000.  OOO, 
9000. OOO, 
9000.  OOO, 

9000. 000, 
203.  OOO, 
9000. OOO, 
-114.  700, 
9000. OOO, 
9000.  OOO, 

9000. OOO, 
1153.  000, 
639.  700, 
9000.  OOO, 
9000. OOO, 
0.  OOO, 

176.  700, 
-311. OOO, 
9000.  OOO, 
9000. 000, 
9000.  OOO, 
9000.  OOO, 

132.  100, 
-262. 600, 
9000. 000, 
9000.  OOO, 
9000. OOO, 
9000. 000, 

42. 730, 
1.  110, 
9000. OOO, 
9000. OOO, 
9000. OOO, 
9000. COO, 

60. 820, 
9000. OOO, 


9000. OOO, 

247. 300, 

523. 000, 

29.  860, 
9000. 000, 
-246. OOO, 
9000. OOO, 

-50. 710, 
9000. OOO, 
9000. 000, 
9000.  OOO, 
9000.  OOO, 
9000. OOO, 

-225. 300, 
-222. 700, 
9000. OOO, 
9000. OOO, 
9000. OOO, 
9000. OOO, 

9000. OOO, 
9000. OOO, 
9000. OOO, 
9000. OOO, 
9000. OOO, 
9000. OOO, 

-80. 480, 
9000.  OOO, 
9000. OOO, 
9000. OOO, 
9000. OOO, 
0.  OOO, 

-17. 780, 
9000. OOO, 
9000. OOO, 
9000. OOO, 
9000. OOO, 
9000. OOO, 

59.  160, 
9000. OOO, 
9000. 000, 
9000. OOO, 
9000. OOO, 
9000. OOO, 

29. 770, 
9000. OOO, 


9000. OOO, 

453. 400, 
82. 770, 

9000. 000, 
-241. 700, 
0. OOO, 
9000. OOO, 

-30. 280, 
9000. OOO, 
9000. OOO, 
9000. OOO, 
9000. OOO, 
9000. OOO, 

-124. 600, 
9000. 000, 
9000. OOO, 
-906. 500, 
9000. OOO, 
9000. OOO, 

9000. OOO, 

417. 400, 
9000. OOO, 

9.  OOO, 
9000. OOO, 
9000. OOO, 

9000. OOO, 
9000. OOO, 
9000. OOO, 
9000. OOO, 
9000. OOO, 
9000. OOO, 

9000. OOO, 
9000. OOO, 
9000. OOO, 
9000. OOO, 
9000. 000, 
0. OOO, 

26. 590, 
9000. OOO, 
9000. OOO, 
9000. OOO, 
9000. OOO, 
9000. OOO, 

55.  970, 
9000. OOO, 


-289.  300 


-99. 560 


-319. 200 


837. 900 


9000. OOO 


9000. OOO 


9000. OOO 


9000. OOO 


*  9000.000,  9000.000,  9000.000,  9000.000, 

*  9000.000,  9000.000,  9000.000,  9000.000, 

*  9000.000,  9000.000,  9000.000,  9000.000, 

*  9000.000,  9000.000,  9000.000,  9000.000, 

*  9000.000,  0.000/ 

DO  5  1  =  1, 32 

AI ( I,  i ) =  AI < I  ) 

A I  <  I,  2) =A2<I ) 

AI < I,  3 ) -A3  <  I > 

AI  <  I#  4 }  =  A4  ( I  ) 

AI (I, 5) =A5 ( I > 

AI < I, 6) =A6  C I  ) 

AI  ( 1,  7 )  =  A7  <  I ) 

A I < 1 ,  8 )  = AS ( X ) 

AI (I,  S)*A9<X> 

AKI,  10)=A10<I) 

AI  (I,  11 )  =Ali  <  I ) 

AKI,  12  )  =  A12  <  I  ) 

Aid,  13  )  =  A13  ( I ) 

Aid,  14 )  «  A14  CD 
AI (I,  15 ) - A 1 5  < I ) 

AI il,  IQ) =A16(I ) 

AI (I,  17 ) -A17(I ) 

AKI,  18  )  SA18  <  I ) 

AKI,  19 )  =  A19  <  I ) 

Aid,  20  )  *  A20  ( I ) 

AKI,  21 )  =A2l  ( I ) 

AKI,  22)=A22d) 

AI ( I,  £3 ) =A23 ( I ) 

AI <1,  24) =A24< I ) 

AKI,  25) =A25<  I  ) 

AKI,  26  )  =A26  <  I ) 

AI Cl,  27 ) =A27  C I ) 

AKI,  28  )  =A28  (  I  > 

AI C I,  29) =A29<I ) 

AI ( I,  30) =A30CI ) 

AI Cl,  31 ) =A31 < I ) 

AI C I,  32) =A32<I ) 

5  CONTINUE 

DO  10  1=1, 57 
SKI)  =  RR  ( I  ) 

OKI)  =QQd  ) 

MAIN ( I ) =MAINSG< I ) 

10  CONTINUE 
RETURN 
END 


9000. 000, 
9000. 000, 
9000. 000, 
9000. 000, 
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oon  ono  o  n  ona 


C 

c  . . . . . . . . . . . . 

c  * 

c  #  <<<  MAIN  PROGRAM  TO  CALL  UNIFAC  SUBROUTINES  >>> 

C  * 

C  * 

C  *  RTI  INTERACTIVE  VERSION!  10-8-86  (WRITTEN  BY  TONY  ROGERS) 

C  * 

C  . . . . * . . . . . 

c 

c 

IMPLICIT  REAL*3(A~H,  0-2) 

CHARACTER *25  NCOMP 

DIMENSION  X ( 10 ) , ACT ( 10 ) , DACT  < 10, 10 ) ,  TACT  < 10 ) , DLACT ( 10 ) , PVAP ( 10 ) , 
*  AI <  40, 40 ) ,  RI ( 76 ) , 01 ( 76  > , MAIN  <  76 ) ,  NCOMP  < 10 ) ,  HENRY ( 10 ) 

DATA  MDL  /3/ 

DATA  LIQ  /O/ 

DATA  NVAP  /!/ 

DATA  NC  /2/ 

DATA  IPASS  /O/ 

DATA  IFLAG  /O/ 

DATA  NH  /O/ 

-  SET  UP  I/O  FILE  NAMES  - 

OPEN < 6, FILE* ' UNIFAC.  DAT ' ,  STATUS* ' OLD ' ,  ACCESS* ' SEQUENTIAL '  , 

♦FORM* 'FORMATTED' ) 

OPEN <  7, FILE* ' UNIFAC.  OUT ' ,  STATUS* ' NEW ' ,  ACCESS* ' SEQUENTIAL ' , 

♦FORM* ' FORMATTED ' ) 

OPEN  <  8, FILE* ' IOWN. DAT ' , STATUS* ' OLD ' ,  ACCESS* ' SEQUENTIAL ' , 

♦FORM* ' FORMATTED ' ) 

1  READ < 6,*)  MODEL,  IOUT,  NDIF,  NACT,  LIQ 
IF( IPASS. GT. 0)  GOTO  2 

IF ( IOUT. GT. 0 )  OPEN ( IOUT,  FILE* ' BACKUP. OUT ' , STATUS* ' NEW ' , 

♦ACCESS* 'SEQUENTIAL', FORM* ' FORMATTED ' ) 

-  PRINT  "START-UP*  MESSAGES  - 

2  IF<LIQ. GT. 0>  GOTO  5 
WRITE ( 7, 901) 

WRITE ( *, 901) 

GOTO  lO 
5  WRITEC7, 902) 

WRITE (  •,  902) 

-  INTERACTIVE  READ  STATEMENTS  - 

10  WRITE (*,  *<1X,A\)')  'Enter  database  ID;  C1J  for  VLE,  C23  for  LLE, 
&r  C3J  for  "Environmental"  ==>  ' 

READ ( *, ' (BN, 13) ' )  MDL 

IF( (MDLrNE. 1 ). AND. (MDL. NE. 2) . AND. (MDL. NE. 3) )  GOTO  10 
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nno 


WRITE ( *,  '</>') 

C 

IF<  LIQ. GT. 0 )  GOTO  100 
C 

20  WRITEC *,'( IX, AN  )' )  'Enter  total  number  oi  components  =  =  >  ' 

READ  (•*,  '  (BN,  13)  '  )  NC 
IFCCNC.  LE.  0).QR.  CNC.  GT.  10))  GOTO  20 
WRITE C *, '(/)') 

IF  <  HC. GE. 2 )  GOTO  30 

-  ERROR  MESSAGE  FOR  HUMBER  OF  COMPONEHTS  - 

WRITEC*,  'CSX,  A)')  '*  Number  of  components  must  be  2  or  greater  1  *' 

WRITEC  *,  '(//>') 

GOTO  130 
C 

30  WRITEC *,'( IX, AN  )' )  'Enter  absolute  temperature  In  CK3  ==>  ' 

READC  *,  *  >  T 

IF  <  T.  LE.  O.  DO  )  GOTO  30 

WRITEC  *,  '(/)') 

C 

DO  35  L=i, HC 

WRITEC  *,  '  C  SX,  A,  12,  AN  )  '  >  'Name  C <25  chars.)  for  Component ', L, 

& '  s=>  » 

READC *, ' CBN, A) ' )  NCOMPCL) 

WRITEC *,  * C 5X,  A,  12,  A\ ) ' )  'Liquid  mole  fraction  for  Component ', L, 

1'  »*>  ' 

READC*, *)  XCL) 

WRITEC*,  '(/)') 

35  CONTINUE 

WRITEC  *,  '</)') 

C 

40  WRITEC *,' C IX, AN )' )  'Calculate  K-factors  Cy/x)  for  components?  Cl  =  \ 
&,  0=N  3  =  *>  ' 

READC  *,  ' CBN, 13) ' )  NH 

IFC CNH. NE. 0) . AND.  CNH. NE. 1 )  )  GOTO  40 

WRITEC*,  '  U)  '  ) 

IF CNH. EQ.O)  GOTO  100 
C 

50  WRITEC C IX, AN )' )  'Enter  [13  far  vapor  pressure  entry,  C23  for  Ar 
ttoine  coefficients  =*>  ' 

READC*, ' CBN, IC> ' )  NVAP 

IFC CNVAP. NE. 1) . AND. CNVAP. NF- 2) )  GOTO  50 
WRITEC  *,  ' C / )  '  ) 

IF CNVAP. EQ. 2)  GOTO  60 
C 

DO  55  K=l, NC 

WRITE  C  *,  ' <  5X, A, 12, AN  >  ' )  'Enter  Component ' , K,  '  v , p.  in  mm  Hg  » *>  ' 
READ C  *,  * )  PVAPCK) 

55  CONTINUE 
C 

WRITEC  *,  '  C  /  )  '  ) 

GOTO  100 
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non  n  noon  non 


60  DO  70  K=i,NC 

WRITE<*,  '  (5X,  A,  12,  A\ )  '  )  'Antoine  Coefficient  A  (v.p.  in  mm  Hg>  for 
&  Comp.  ' ,  K,  '  =  =>  ' 

READ (  *,  *  )  ANT  A 

WRITER  '  <  SX,  A,  12,  A\ )  '  )  'Antoine  Coefficient  B  (v.p.  in  mm  Hg  )  for 
&  Comp.  ' ,  K,  '  *  =  >  ' 

READ(  *, * )  ANTB 

WRITE  <  *,  '  <  5X,  A,  12,  A\ )  '  )  'Antoine  Coefficient  C  (v.p.  in  mm  Hg  )  for 
L  Comp. ' , K, '  ==>  ' 

READ (  *,  *  )  ANTC 

PVAP <  K ) =D£XP <  ANTA- ANTB/ ( T  + ANTC > ) /760. DO 
WRITE  <  <*,  '(/)') 

70  CONTINUE 

100  IF ( MDL. EQ. 1 )  CALL  UVLE( RI, QI,  AI, MAIN ) 

IF ( MDL. EQ. 2 >  CALL  ULLE ( RI, QI, AI, MA IN ) 

IF ( MDL. EQ. 3 )  CALL  ENVIRONC RI, QI, AI, MAIN ) 

IF ( LIQ. GT. 0 )  GOTO  120 

CALL  FINOUT ( NC, RI,  QI, AI, MAIN, NG > 

REWIND  6 

WRITE  <  *,  '<//)') 

PAUSE 

WRITE<  *,  '(//)') 

CALL  PARAM<NC, NG, T) 

CALL  UNIFA  <  NDIF,  NACT,  NC,  NG,  T,  X,  ACT,  DACT,  TACT ) 

-  CALCULATE  K-VALUES  CY/X3  FOR  ALL  COMPONENTS  (NH=1>  - 

WRITE  <  7, 903 )  T 
WRITE <  *, 903 )  T 
DO  110  K*l, NC 

HENRY ( K ) *-999. DO 

IF ( NH. GT. O )  HENRY  <  K ) = ACT ( K ) *PVAP ( K ) 

WRITE ( 7, 904 )  K,  NCOMP ( K ) , X( K) , ACT ( K ) , HENRY ( K ) 

WRITE ( ♦, 904)  K,  NCOMP < K ),  X ( K ) ,  ACT <  K ) ,  HENRY ( K ) 

WRITE ( 7, 905) 

WRITE ( *, 905) 

110  CONTINUE 

GOTO  ^30 

-  CALCULATION  OF  AQUEOUS  SOLUBILITIES  (LIQ=i)  - 

120  CALL  XLIQUIDt NC, MDL, MODEL,  lOUT,  NDIr,  NACT,  RI,  QI, AI, MAIN) 

130  CONTINUE 

-  USER  PROMPT  TO  RE-START  PROGRAM  OR  END  - 

WRITE <*,'( //IX,  A\ )' )  'Enter  £1]  to  re-start  program  or  tQ3  to  end 
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non 


&==>  ' 

READ< *,  ' < BN, 13)  '  >  IFLAG 

IF< < IFLAG. NE. O) . AND. < IFLAG. NE. 1 > >  GOTO  130 
WRITE ( *,  '  (/////  )  ' ) 

IF (IFLAG. EQ. 0)  GOTO  999 

REWIND  6 

IFLAG=Q 

IPASS=i 

GOTO  1 

999  CONTINUE 
C 

STOP 


---  FORMAT  STATEMENTS  - 


C 

C 

C 

C 

C 

C 

C 

C 

c 

c 


c 


c 


c 


901  FORMAT < 20X, ' ***  ACTIVITY  COEFFICIENT  CALCULATION  ***',///, 

&10X, 'NOTE:  Vapor  pressure  data  (or  Antoine  coefficients)  must  be  s 
iupplied',  /,  16X,  'for  all  components  to  calculate  their  Henry's  Con 
istants ' , / // ) 

902  FORMAT < 15X, ' ***  BINARY  LIQUID-LIQUID  EQUILIBRIUM  ROUTINE  ***',///, 
&10X, 'NOTE:  Enter  the  most  dilute  chemical  as  Component  #1. ',///> 

903  FORMAT (25X. 'TEMP.  ■ ' ,  F8.  2,  '  CK] ' , //,  5X,  ' ID  #',  5X,  ' COMPONENT  NAME', 
&4X,  'LIC.  M-FRAC',  5X,  'ACT.  COEFF.  ' ,  5X,  '  H  <TU>',//) 

904  FORMAT  ( bX,  12,  6X,  A1S,  2X,  E13.  5,  3X,  F12.  5,  3X,  F8.  1 ) 

905  FORMAT <///) 

END 


<<<  SUBROUTINE  TO  HANDLE  INPUT/OUTPUT  OF  DATA  >>> 


* 

* 


SUBROUTINE  FINOUT(NC, RI, QI, AI, MAIN, NG) 


IMPLICIT  REAL *8  <  A-H, 0-2  > 

COMMOM/UNIF/RT ( 10,  10),  QT(10,  10),  TALK  10,  10 ) ,  S ( 10,  10 ) ,  F ( 10 ) , Q ( 10 ) , 

*  R  < 10 ) , P ( 10, 10  > 

DIMENSION  AI ( 40,  40 ) ,  RI ( 76 ) ,  QI < 76 ) ,  MAIN ( 76 ) , 

*  NGM< 10), MS< 10,  10, 2), NY (10, 20 ) , JH ( 76 ) , IH ( 20 ) 

IF( IOUT. EQ.  O)  I0UT*7 
NK  =  NC 


DO  10  1=1, 10 
DO  10  J=l, NK 

QT (  I,  J  )  =0.  DC 
RT ( I,  J) =0. DO 
10  CONTINUE 

IF( MODEL. NE. 1 )  GO  TO  30 
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c 


NG=NK 


DO  20  1*1,  NK 

READ  (8,  *)  P.T(  I,  I  >,  QT<  I,  I  >,  <  P  ( I,  J ) ,  J*l,  NK ) 

20  CONTINUE 
30  CONTINUE 
C 

IF < MODEL. EO.  1)  GO  TO  290 
READ  CS,  * )  IOWNRQ, IOWNP 
IF ( IOWNRQ.  EQ.  0  >  GO  TO  50 
C 

DO  40  1*1, IOWNRQ 

READ  <  8,  * )  K,  RI  <K>,  QI  <K) 

40  CONTINUE 
C 

50  IF < IOWNP. EQ. 0)  GO  TO  70 

✓"* 

DO  60  1*1, IOWNP 

READ  <  8,*)  J,  K,  AI<J,K> 

60  CONTINUE 
70  CONTINUE 
C 

DO  80  1*1, NK 
DO  30  J*l, 10 
DO  80  K*l, 2 
MS(I, J,K)»0 
30  CONTINUE 
C 

DO  90  1*1,  76 
JH  < I ) *0 
90  CONTINUE 
C 

NMAX*0 

DO  100  1*1, NK 

WRITE < *, ' < IX, A, 12, A\ ) ' )  'Enter  number  of  types  of  subGroups  in 
(component ' , I, '  **>  ' 

READ  <  *,  ' ( BN, 13 ) ' )  NGRP 
WRITE <  *,  '</)') 

DO  95  KR* 1, NGRP 

WRITE (*, '(5X,A,I2,A\)')  'Enter  ID#  of  aubGroup ' , KRf '  **>  ' 
READ<  '  <BN,  13)  '  >  MS(I,KR,1) 

WRITE<  *,  '  (  5X,  A,  13,  A,  12,  A\ )  '  )  'Enter  number  of  aubGroup ',  KR, 

'  in  component ', I, '  *=>  ' 

READ  (  *,  '  <  BN,  13  >  '  )  MS<I,KR,2> 

WRITE ( *,  '</)') 

CONTINUE 
WRITE<*, '(/)') 

IF ( NGRP. GT. NMAX )  NMAX=NGRP 
100  CONTINUE 
C 


IC*1 
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c 

DO  ISO  1=1, NK 
C 

DO  150  J=l, NMAX 

IF<MS( I, J, 1). EQ. O)  GO  TO  ISO 
'  IH ( IC ) =  MS (I, J, 1) 

IF  < IC. EQ. 1)  GO  TO  140 
ZP(IH(IC).EQ.  IH(IC-D)  00  TO  150 
IPdHdO.GT.  IH<IC-D)  GO  TO  140 
IF< IC. GT. 2)  GO  TO  110 
IHH=IH< 1 ) 

IH< 1 >  =  IH<  2) 

IH ( 2 ) =  IHH 
GO  TO  140 
110  I1=IC-1 

C 

DO  130  12  =  1,11 

IFC IH< IC ) . GT. IH  < 12) )  GO  TO  130 
IF< IH< IC  > . EQ. IH  < 12  >  >  GO  TO  150 
I4=IC- 12 
C 

DO  120  13=1,14 

IH ( IC+I-I3 )=IH<IC-I3> 

120  CONTINUE 

C 

IH< 12) =  MS( I,  J,  1) 

130  CONTINUE 

C 

140  IC*IC+1 

IF(IC. GT. 20)  WRITE <7, 905) 

IF < IC. GT. 20 )  WRITE <  *, 905 ) 

IF  < IOUT . NE. 7. AND. IC. GT. 20)  WRITE < IOUT, 905) 
150  CONTINUE 
C 

160  CONTINUE 
C 

IC-IC-1 

C 

DO  170  1=1, IC 
JH  < IH ( I ) >=I 
170  CONTINUE 
C 

DO  180  1  =  1,10 
DO  180  J=l, 20 
NY ( I, J)=0 
180  CONTINUE 
C 

DO  200  1=1, NK 
C 

DO  190  J=i. 10 

IF  <  MS  < I, J, 1 ) . EQ. 0 )  GO  TO  200 
N1=MS< I, J, 1) 

N2=MS( I, J, 2) 


80 


IFCNi.EQ.  0)  GO  TO  200 
N3=JH(N1 ) 

NY  <  I,  N3  )  =  N2 
190  CONTINUE 
C 

200  CONTINUE 
C 

1*0 

NGMGL*0 

C 

DO  210  K=l, IC 
NSG*IH(K> 

NGJINY  =  MAIN  ( MSG ) 

IF  ( NGMHY.  NE.  NGflGL )  1  =  1  +  1 
NGM  < I >  *NGMNY 
NGMGL*NGMNY 
DO  210  J=1,NK 

RT ( I, J ) =RT ( I, J) +NY< J, X ) *RI ( NSG  > 
QT(I, J)=GT<I, J)+NY<J, K) »QI(NSG) 

210  CONTINUE 
C 

NG=I 

WRITE ( 7, 906 )  < IH < K ) , K=l, IC ) 

WRITE ( *, 906 )  < IH ( K ) , K=i, IC ) 

WRITE <7, 907)  <MAIN<IH<K) ),  K*l, IC) 
WRITE (  *,  907)  < MAIN <  IH ( K )  >,  K»l,  IC) 
WRITE <7, 90S) 

WRITE <  *,  90S) 

C 

DO  220  1*1, NK 

WRITE  (  *,  909 )  I,  <NY(I,  K),  K=l,  IC) 
WRITE(7,909)  I,  < NY  <  I,  K) ,  K=l,  IC ) 

220  CONTINUE 
C 

WRITE (7, 913) 

WRITE ( +  , 913) 

IF  < IOUT. EQ. 7 )  GO  TO  240 
WRITE (IOUT, 906)  < IH(K ) , K=l, IC ) 

WRITE ( IOUT, 907)  <HAIN<IH(K) ),  K  =  l, IC) 
WRITE (IOUT, 908) 

C 

DO  230  1=1, NK 

WRITE< IOUT, 909)  I,  ( NY ( I,  K ) , K  =  1, IC) 
230  CONTINUE 
C 

WRITE( IOUT, 913) 

240  CONTINUE 
C 

DO  250  1=1, NG 
DO  250  J=l, NG 
NI=NGM(I) 

N J  =  NGM  <  J  ) 

P(I,  J  )  *AI  <  NJ,  NI) 
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c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


250  CONTINUE 

WRITE ( 7, 910) 

WRITEf  *, 910) 

DO  260  K*  1, IC 
NN- IH ( K ) 

WRITE  (  *,  911 )  NN,  RI ( NN ) , QI ( NN  > 
WRITE ( 7, 911 >  NN,  RI <  NN ) ,  QI <  NN ) 

260  CONTINUE 

WRITE <  7, 913) 

WRITER  *, 913) 

IF ( IOUT. EQ. 7  >  GO  TO  230 
WRITE< I0UT, 910) 

DO  270  K  =  i, IC 
NN-IH(K) 

WRITE; IOUT, 911)  NN,  RI(NN), QI(KN) 
270  CONTINUE 

WRITE; IOUT, 913) 

230  CONTINUE 
290  CONTINUE 

WRITE (7, 902) 

WRITE; *,902) 

DO  300  1=1, NG 

write;*, 901)  ; p ; i,  j > ,  j=i,  ng > 
write;7,9oi>  ;p;i,  j),  j=i,  ng) 

300  CONTINUE 

WRITE ( 7, 913) 

WRITE; *, 913) 

IF ; MODEL. EQ. 0 )  WRITE ; 7,  903  > 

IF ; MODEL. EQ.  0)  WRITE; *,903) 

IF ; MODEL. EQ. 1)  WRITE; 7, 912) 

IF< MODEL. EQ.  1 )  WRITE;*, 912) 

IF (IOUT. EQ. 7)  GO  TO  320 
WRITE (IOUT, 902) 

DO  310  1=1, NG 

WRITE; IOUT, 901 )  <P< I, J), J=i, NG) 
310  CONTINUE 

WRITE; IOUT,  913) 

IF( MODEL. EQ. 0)  WRITE < IOUT, S03 ) 

IF< MODEL. EQ. 1 )  WRITE ( IOUT, 91 2 ) 

320  CONTINUE 


DO  330  1=1, NK 


ft 
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Q( I > *0. DO 
R  < I ) -0. DO 
DO  330  K=l, KG 

Q ( I ) *Q ( I ) ♦QT  <  K, I) 

R(I)*R<I>-RT<K, I) 

330  CONTINUE 

U 

DO  340  1*1, NK 

WRITE<  *, 904 )  I,R(I),Q(I) 

WRITER  7, 904 )  I,R(I),Q(I) 

340  CONTINUE 
C 

IF ( IOUT. EQ.  7 )  GO  TO  360 
DO  350  1*1, NK 

WRITE ( IOUT, 904)  I,R<I),Q(I> 

350  CONTINUE 
360  CONTINUE 

-  FORMAT  STATEMENTS  - 

901  FORMAT  < IX, 1QF12. 3 ) 

902  FORMAT < 2X, ' INTERACTION  PARAMETERS',/) 

903  FORMAT < IX, 'UNIFAC  MOLECULAR  R  AND  Q',/) 

904  FORMAT  ( 15, 2F1S. A ) 

905  FORMAT < IX,  * **  WARNING:  NO,  OF  SUBGROUPS  MUST  NOT  EXCEED  20  **') 

906  FORMAT </,  IX,  'SUB  GROUPS  :',20I3) 

907  FORMAT < IX,  'MAIN  GROUPS 2013 ) 

908  FORMAT < IX, 'COMPONENT' ) 

909  FORMAT (6X, 12, 5X, 2013) 

910  FORMAT (IX, 'GROUP  R-  AND  Q-VALUES',/) 

911  FORMAT< IX, 13,  2F10. 4) 

912  FORMAT < IX, 'SPECIFIED  UNIOUAC  R  AND  Q',/) 

913  FORMAT ( // ) 

RETURN 
END 

C 

C 

. . . . . 

c 

C  <<<  UNIFAC  BINARY  LIQUID-LIQUID  FLASH  ROUTINE  >>> 

C 

C 

C  NEWTON  -RAPHSON  ALGORITHM  ( GAUSS- JORDAN  MAXIMUM  PIVOT  STRATEGY) 

- - 

C  THIS  IS  A  SUBROUTINE  TO  IMPLEMENT  THE  NEWTON-RAPHSON  ALGORITHM  FOR 

C  SOLVING  SYSTEMS  OF  NONLINEAR  ALGEBRAIC  EQUATIONS.  A  VARIATION  OF 

C  THE  GAUSS-JORDAN  MAXIMUM  PIVOT  STRATEGY  IS  EMPLOYED  TO  DETERMINE 
C  THE  INVERSE  OF  THE  JACOBIAN  MATRIX.  THE  CORRECTION  FACTORS  ARE 
C  CALCULATED  IN  AN  ITERATIVE  MANNER  TO  BRING  THE  ADJUSTABLE  VARIABLES 
C  WITHIN  A  SPECIFIED  TOLERANCE. 

C 

C  M  *  NUMBER  OF  COLUMNS  IN  MATRIX  C 
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c  H  =  NUMBER  OF  ROWS  IN  MATRIX  C 

C 

C 

. . . . . . . . . . 

SUBROUTINE  XLIQUID(NC,  MDL,  MODEL,  IOUT,  NDIF,  NACT, RI, QI, AI, MAIN) 

C***# ********************************************************  *********** 

IMPLICIT  REAL*8< A-H, O-Z) 

CHARACTER *25  NAME 

DIMENSION  C  < 10, 11 ) , KNEW < 10 ) ,  DX < 10 > ,  IR < 10 ) ,  F ( 10 ) ,  XMW < 10 ) , 

XI < 10  >  ,  X2< 10 )  ,  XOLD < 10  ) ,  FOLD ( 10 ) ,  TEST  < 10 ) ,  IE ( 10  > , 

ACT1 ( 10) , DACT1 < 10,  10 ) ,  TACT! < 10  > ,  DLACT1 <  10 ) , 

ACT2< 10), DACT2< 10,  10  >  ,  TACT2 ( 10  > ,  DLACT2  < 10), 

XE < 10, 2), TEMP (10),  FMINC10,  2),  NAME(IO), 

AI ( 40, 40  > ,  RI <  76 ) ,  QI ( 76 ) , MAIN  <  76 ) 

C 

C  -  SET  LENGTH  AND  WIDTH  OF  JACOBIAN  MATRIX  - 

C 

N=NC 

M=N*I 

INDEX=0 

C 

C  -  INITIALIZE  ARRAYS,  VECTORS  - 

C 

DO  5  1=1, 10 
C 

DO  6  J*i,N 

XE( I, J ) =0. DO 
FMIN< I, J ) =0. DO 
6  CONTINUE 
C 

TEMP  < I ) =0. DO 
NAME  < I )  =  '  ' 

5  CONTINUE 
C 

C  -  BEGIN  OUTER  LOOP  TO  READ  MULTIPLE  DATA  SETS  - 

C 

READ  (6,*)  IMAX,  TOL 
C 

WRITE <*,'< IX, A\ )' )  'Enter  number  of  data  sets  (binary  pairs)  *=>  ' 
READ ( *, ' <BN, 13) ' >  NDS 
WRITE ( *, '</)') 

C 

DO  777  KN=i, NDS 
C 

C  -  READ  DEFAULT  PARAMETERS  - 

C 

WRITEC *, ' < IX, A, 12, A\ ) ' >  'Give  name  ( <25  chars.)  of  binary  system', 
4XN, '  =»>  ' 

HEAD< *, ' < BN, A) ' )  NAME ( KN ) 

WRITE  < *,'</)') 

C 

WRITE (*,'( IX, A\ )' )  'Enter  absolute  temperature  in  CK3  ==>  ' 

READ ( *,  * )  TEMP ( KN ) 
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T*T£MP  <  KN ) 

WRITE ( *, '</)') 

C 

WRITE  <  7, 902 )  KN, NAME ( KN ) , TEMP ( KN ) 

WRITE*  *, 902)  KN, NAME ( KN ) , TEMP ( KN ) 

WRITE <7, 903) 

WRITE ( *, 903) 

C 

DO  7  LP=i,  N 

WRITE ( *,  ' < 5X,  A,  12,  A\ ) ' >  'Enter  molecular  weight  for  component ' ,  LP, 
&'  ==>  ' 

READ ( *, * )  XMW  <  LP ) 

C 

WRITE  (  *,  '  <  5X,  A,  12,  A\ )  '  )  'Enter  composition  guess  for  component', 
&LP, '  *=>  ' 

READ  <  *, *  >  XNEW  <  LP ) 

WRITE ( *,  '</)') 

7  CONTINUE 

WRITE* *,  '*/>') 

C 

CALL  PINOUT  ( NC,  RI,  QI,  AI,  MAIN,  NG ) 

REWIND  8 

- INITIALIZE  COEFFICIENT  MATRICES  ( C, FOLD, XOLD,  DX,  IE ) - 

DO  10  L»1,H 

FOLD  *  L ) *0. DO 
XOLD  *  L  > 3  XNEW  <  L ) 

DX  <  L ) =0. DO 
IE<L)»0 

DO  20  LL=1, M 

C(L,  LL)=0. DO 
20  CONTINUE 

10  CONTINUE 

-  WRITE  "WAIT"  MESSAGE  TO  SCREEN  - 

WRITE (7, 904) 

WRITE <  *, 904) 


<<<  START  OF  NEWTON-RAPHSON  ITERATION  LOOP  >>> 


DO  939  ITER=1, IMAX 
INDEX3 ITER  - 1 

WRITE<7, 905)  INDEX 
WRITE <  *, 905 )  INDEX 

-  EVALUATE  FUNCTION  VECTOR  <F)  -- 
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C 

XI < I ) =  XNEW  < 1 ) 

XI  <  2>  =  1. DO-XK i) 

X2  < 1 ) =XNEW ( 2 ) 

X2<2>=1. D0-X2( 1) 

C 

CALL  PARAM  <  N, MG, T ) 

CALL  UNIFA  <  NDIF,  NACT,  N,  NG,  T,  XI,  ACT1,  DACT1,  TACT1 ) 

CALL  UNIFA  (  NDIF,  NACT,  N,  NG,  T,  X2,  ACT2,  DACT2,  TACT2  > 

C 

DO  21  J  =  l, N 

F  <  J) =Xi ( J ) *ACT1 ( J ) -X2  ( J  > *ACT2 ( J ) 

TEST <  J )  =  DSORT  < ( FOLD < J ) )**2+<F< J) )**2> 

21  CONTINUE 

-  TEST  FOR  CONVERGENCE  OF  SOLUTION  - 

LQGIC=-i 

DIFF=DSQRT( ( KNEW < 1 > -XOLD ( 1 > >  *  *2-*- <  XNEW  (  2 ) -XOLD  ( 2  > )**2> 
IF ( DIFF. GT. TQL )  LOGIC=l 

DO  22  1=1, N 

IF (TEST (I ). GT. TOL)  LQGIC=i 
IF<DX(I). GT. TOL)  LOGIC* I 

22  CONTINUE 

IF ( LOGIC. LT. 0)  GOTO  d 

DO  25  1*1,  N 

FOLD  < I ) =F  < I ) 

XQLD  < I ) =XNEW  < I ) 

25  CONTINUE 

-  CALCULATE  PARTIAL  DERIVATIVES  IN  JACOBIAN  - 

C(  1,  I  )  *  ACT1  <  1  )  -*-XNEW  <  1 )  *DACTi  ( 1 ,  1) 

C ( 1 ,  2 ) = - ACT2  < 1 ) -XNEW ( 2 ) *DACT2 (1,1) 

C(2,  1 ) =-ACTl <  2 ) * ( 1.  DO-XNEW ( 1 ) )*DACT1(2, 1) 

C ( 2,  2 ) *  ACT2 ( 2 > - ( 1. DQ-XNEWC2) )*DACT2(2,  1) 

-  FINISH  LOADING  C  MATRIX  WITH  F  VECTOR  - 

DO  30  1  =  1,  N 

C(I, M)*-F(I) 

30  CONTINUE 

C 

C  **•  GAUSS-JORDAN  ALGORITHM  *** 

C 

C 

C  -  INITIALIZE  ALL  VECTORS  AND  MATRICES  - 

C 

DO  40  1  =  1,  N 
IR  < I ) =0 


86 


oonoon  n  o  nono  nnnnn  n  non 


40 


C 


DX  Cl)  *0.  DO 

JJ=G 

JM=0 

CONTINUE 

DO  50  K=i,N 
PK=0. DO 

-  LOCATE  PIVOT  ELEMENT  - 

DO  60  1=1, N 

IF (I. EQ. IR  <  I ) )  GOTO  60 


70 


DO  70  IK=1,N 

P=DABS( C ( I,  IK) ) 

IF <  P. LT.  PK )  GOTO  70 
PK=P 
JJ*I 
JH=IK 
CONTINUE 


60  CONTINUE 

IR  <  JJ  >  3  J J 

---  NORMALIZATION  STEP  - 

DO  60  JR=1,M  • 

IF<JM. EQ. JR)  GOTO  8G 
C  <  J  J,  JR ) =C  <  J J,  JR )  /C  ( J J,  JM ) 
80  CONTINUE 

C(  JJ,  JM)  =1.  DO 

-  REDUCTION  STEP  - 

DO  90  1=1, N 

IF< I. EQ. JJ)  GOTO  90 


DO  100  JR=1,  M 

IF(JR.EQ.  JM)  GOTO  100 
C  <  I,  JR  >  »C  <  I,  JR  )  -C  <  I,  JM  >  *C  <  J  J,  JR  ) 
100  CONTINUE 


C<I, JM)=0. DO 
90  CONTINUE 

50  CONTINUE 

***  END  GAUSS-JORDAN  MAXIMUM  PIVOT  ROUTINE  *** 


-  RECOVER  SOLUTION  VECTOR  -- 
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c 


DQ  110  1  =  1,  N 


C 

DO  120  J=l,  N 

IF(  (C<I,  J).  LT.  1.  D0>.  OR.  <C(I,  J).  GT.  1.  D0>  )  GOTO  120 
DX< J) =C<  Z, K) 

120  CONTINUE 

C 

110  CONTINUE 

-  CORRECT  ELEMENTS  OF  X  VECTOR  - 

DO  130  1  =  1,  N 

XNEW { I ) =XNEW  < I ) ♦DX  < I ) 

IF ( XNEW  < I ) . LT. 0. DO )  XNEW ( I >  =0. DO 
IF(XNEW(I>.  GT.  1.  DO)  XK£W<I)=1.D0 
130  CONTINUE 

999  CONTINUE 

GOTO  9 

«<  END  OF  ITERATION  LOOP  »> 


-  PRINT  FINAL  MATRIX  <C>  AND  SOLUTION  VECTOR  - 


8  WRITE ( *, 906 )  INDEX 
WRITE <7, 906)  INDEX 


DQ  140  I"  l,  N 

WRITE <,  /,  907)  <C(  I,  J),  J=l,  N>,  C<  I,  M> 

WRITE (  *,  907)  <C(I,J),J=i,N),C(I,M> 
140  CONTINUE 


WRITE  <7, 908 ) 

WRITE  <  *, 908 ) 

DO  150  1  =  1,  N 

WRITE ( 7, 909 )  I ,  XNEW ( I ) , I, F ( I ) 
WRITE  (  *,  909  )  I,  XNEW  <  I  ) ,  I,  F  <  I ) 
150  CONTINUE 


WRITE (7, 910) 
WRITE <  *, 910) 


-  LOAD  SUMMARY  VECTORS  FOR  FINAL  RESULTS  TABLE  - 


DO  160  1  =  1,  N 

DO  155  J  =  l,  N 

IF  <  J.  NE.  I  )  GOTO  152 
IE< J) *IE< J) *1 

152  IF  <  XNEW  <  J ) . GT. XNEW  < I )  )  IE ( J ) =IE( J ) ♦! 


88 


non  o  n  n  o  o  non  ooonono  nnn 


155  CONTINUE 

FMIN  < KM, I >  =F  < I ) 

ISO  CONTINUE 

SORT  COMPOSITION  SUMMARY  VECTOR  - 

DO  165  1*1,  N 

XE ( KN, I ) =  XN£W ( IE(  I ) ) 

165  CONTINUE 

-  CONVERT  AQUEOUS  CONCENTRATION  TO  PPMW  UNITS  - 

XE<  KN,  1 >=1000000.  DO/ <1.  DO* ( < 1.  DO/XE<KN,  1>  > -1. DO ) *XMW < 2 > /XMW < 1 >  > 
GOTO  777 

-  CONVERGENCE  FAILURE  MESSAGE  - 

9  WRITE  <  *, 911 )  IMAX 
WRITE  <  7, 911 )  IMAX 
WRITE  <  7, 910) 

WRITE <  *, 910) 

-  LOAD  SUMMARY  VECTORS  WITH  DEFAULT  VALUES  - 

DO  170  1*1,  N 

XE<  KN, I >  =  -999.  DO 
FMINCKN, I >=-999. DO 

IF(XNEW( 1 > . LT. < 1. DO-TQL> >  GOTO  170 
XE(KN, I ) =XNEW  < I ) 

FMIN ( KN, I ) *F < I ) 

170  CONTINUE 

777  CONTINUE 

-  PRINT  TABLE  OF  FINAL  MISCIBILITY  VALUES  - 

WRITE  <  7, 912) 

WRITE ( *, 912) 

WRITE  <  7, 912) 

WRITE  < *, 913) 

DO  ISO  1*1, NDS 

WRITE<7,  914)  NAME  < I > , TEMP  < I > ,  i XE< I, J> , J-l, N > 

WRITEt  *,914)  NAME  <  I  ) ,  TEMP  ( I ) ,  <  XE ( I,  J  > ,  J  =  1,  N  > 

i  ofi  i-nuTTunr 

-  FORMAT  STATEMENTS  - 

901  FORMAT (A25) 

902  FORMAT <//,  5X,  'DATA  SET',  13,  '  :  ',A25,  '  AT',F7.  2,  '  K') 

903  FORMAT  (5X,  Q  <'-'),  6X,  37  <'-'),// > 
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904  F0RMAT(//,  20X,  'PLEASE  WAIT...  CALCULATIONS  PROCEEDING...",/) 

905  FORMAT (30X, 'ITERATION  ',13) 

906  FORMAT(/,  15X,  'FINAL  SOLUTION  MATRIX:  (AFTER',  13,  '  ITERATIONS)', 

&/  ) 

907  FORMAT < 15X, 7E15. 5 ) 

908  FORMAT (/, 5X, 'SOLUTION  VECTOR 15X, ' FUNCTION  VECTOR:',/) 

909  FORMAT <5X,  'XNEW( ',  12,  '  >  =  ',  Eli.  5,  lOX,  'F<  ',  12,  '  )  =  ',E11.5) 

910  FORMAT ( /// ) 

911  FORMAT(/, 15X, 'ALGORITHM  DID  NOT  CONVERGE  AFTER ' , 14, '  ITERATIONS'/ 

912  FORMAT  (/,8X,  'SYSTEM',  7X,  'TEMP  <K)',5X,  'XI  <PPMW,  I)',5X, 

&'X1  (M-F,  II)') 

913  FORMAT (SX, 6  ('-'),  7X,  8  ('-'),  5X,  12  ('-'),  5X,  12  ('-'),/ ) 

914  FORMAT  ( IX,  A15,  5X,  F8.  2,  4X,  F8.  1,  5X,  1E17.  5) 

RETURN 

END 

C 

C 

C  * . . . . . . . . . * . 

c  * 

C  • 

C  *  «<  SUBROUTINE  TO  CALCULATE  THE  ACTIVITY  COEFFICIENTS  >>> 

C  * 

c  * 

c  *****#****•♦••******•*##****#**•*#*****#*#*•#•*■**■»**#**<  »#***#**»»## 

SUBROUTINE  UNIFA  <  NDIF,  NACT,  NC,  NG,  T,  X,  ACT, DACT,  TACT) 


IMPLICIT  REAL*8(A-H, 0-Z) 

COMMON/UNIF/RT ( 10,  10),  QT(10,  10), TAU(10,  10  > ,  SdO,  10 > ,  F  10),  0(10), 

*  R( 10) , P ( 10, 10) 

DIMENSION  X ( 10 ) , GAM ( 10 ) ,  ACT ( 10 ) ,  DACT ( 10,  10 ) ,  THETA  < 10 )  ’HI (10), 
*RI( 10). 01(10),  OIL (10),  RIL(IO),  QID(iO),  ETAL ( 10 ) ,  TACT ( 1  , UC10, 10), 

* V ( 10, 10 ) ,  DETA ( 10 ) , DS ( 10,  10 ) ,  ETA ( 10 ) ,  TETAR ( 10 ) ,  H3 ( 10,  1 
C 

THETS*0. DO 

PHS3Q. DO 
C 

DO  10  1*1, NC 

THETA ( I ) =X( I ) *0 ( I ) 

PHI (I)*R(I)*X(I) 

THETS*THETS* THETA ( I ) 

PHS*PHS*PHI ( I ) 

10  CONTINUE 
C 

DO  20  1*1, NC 

THETA ( I ) =THETA < I ) /THETS 
PHI ( I ) *PHI ( I ) /PHS 
RI(I)*R(I)/PHS 
RIL ( I ) *DLOG ( R I ( I )  ) 

01 ( I ) *0 ( I ) /THETS 

QID ( I ) =1. DO-RI ( I ) /QI ( I ) 

OIL ( I ) =DLQG ( QI ( I ) ) 

20  CONTINUE 
C 


90 


DO  30  1  =  1, NC 

XX  =  F  < I ) *Q ( I ) * ( 1.  DQ-QIL< I ) 5 -RI < I ) -RIL< I ) 

XX=XX-5. D0*0 ( I ) * < QID ( I ) f RIL ( I > -QIL< I > ) 

GAM ( I ) =  XX 
30  CONTINUE 
C 

DO  50  1=1, NG 

TETAR  <  D  =0. DO 
ETA (I )*0. DO 
C 

DO  40  J  =  1 ,  NC 

ETA  < I ) =ETA  <  I)  -*-S<  I,  J  MX<  J  > 

TETAR  <  I )  =TETAR  <  I )  +QT  <  I,  J  MX  <  J  > 

40  CONTINUE 
C 

ETAL<IMDLOG<ETA<U  ) 

50  CONTINUE 
C 

DO  70  1=1, NC 
C 

DO  60  J=i, NG 

U< J, I  MS  <  J, I ) /ETA <  J ) 

V  <  J, I) =U  <  J, I ) *TETAR  <  J ) 

GAM  < I ) *GAM  < I ) - V  <  J,  I ) -QT <  J,  I ) *ETAL<  J ) 

60  CONTINUE 
C 

ACT<IMDEXP<GAM<U  ) 

IF < NACT.  EQ.  1 )  ACT<IMACT<IMX<D 

70  CONTINUE 
« 

C 

IF<NDIF. EQ. 0)  GO  TO  160 
IF<  NDIF. EQ. 2 )  GO  TO  110 
C 

DO  90  1=1, NC 
DO  90  J=l, NC 

XX=Q  < I >  *01 <  J  M < 1. DO-5.  DQ*QID< I ) *QID<  J >)  +  (!.  DO-RI < I) ) * 
&  ( 1. DO-RI <  J  > ) 

C 

DO  80  K=l, NG 

XX*XX  +  U <  K,  I)*<V<K,  J) -QT  <  K, J) ) -U < K, J ) *QT < K, I) 

60  CONTINUE 
C 

DACT  < I, J ) =XX 

DACT  ( J,  IMXX 

IF < NACT. EQ. 1 )  GO  TO  90 

DACT  < I, J ) =DACT  < I, J ) *ACT  < I ) 

IF  <  J. EQ. I)  GO  TO  90 
DACT ( J, U  =DACT  <  J,  I) *ACT  <  J i 
90  CONTINUE 
C 

IF<NACT. EQ. 0)  GO  TO  110 
C 


DO  100  1=1, NC 
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c 

c 


c 


c 


c 


c 


c 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 


DO  iOO  J=i, NC 

DACT  < I, J) =ACT ( I >  *  <  DACT ( I, J ) -1. DO) 

IF ( J.  EQ.  I )  DACT < I, J  >  =DACT  < I,  J ) +  DEXP ( GAM  < I ) ) 

IOO  CONTINUE 
110  CONTINUE 

IF ( NDIF. EQ. 1 )  GO  TO  160 

DO  130  K=i, NG 

DETA  <  K ) =0. DO 
DO  130  1*1, NC 

DS(K, I) *0. DO 

DO  120  M*1,NG 

IF<QT<K,  I).EQ.  0.  DO)  SO  TO  120 

DS<K, I ) -DS(K, I) -QT(M,  I ) *DLQG<  TAU (H,  X ) ) »TAU  ( M,  K ) /T 
120  CONTINUE 

DETA ( K  >  =  DETA <  K ) *DS ( K,  I ) *X  < I ) 

130  CONTINUE 

DO  150  1*1, NC 

TACT  < I ) *0. DO 

DO  140  K=l, NG 

H3<  X, I ) *  < -S<K, I > *DETA<K) /ETA<K >  +DS  ( K,  I > ) /ETA(K ) 

HH*H3<K,  I ) • ( TETAR ( K ) -QT ( K,  I > *ETA < K > /S< K,  I) ) 

TACT  < I ) *TACT ( I ) -HH 
140  <  CONTINUE 

TACT  < I ) *TACT ( I >  * ACT  < I ) 

150  CONTINUE 
160  CONTINUE 
RETURN 
END 


PARAW  CALCULATES  SOME  COMPOSITION- INDEPENDENT  QUANTITIES: 
TAU,  S,  AND  F,  TO  BE  USED  IN  UNIFA.  PARAM  MUST  BE  CALLED 
PRIOR  TO  UNIFA. 


SUBROUTINE  PARAM  <  NC, NG, T ) 

IMPLICIT  REAL*6< A-H, 0-2) 

COMMON/ UNI F/RT (10,  10 ) ,  QT < 10, 10 ) ,  TAU < 10,  10 ) ,  S ( 10,  10 ) ,  F( 10 > ,  Q ( 10 ) , 
R<  10),  P<  10,  10) 


DO  10  1=1, NG 


■HH 
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nnnnnnonnnnnnoonnnnnnnon  ooooonnno 


c 


c 


c 


DO  10  J=1,NG 

TAU ( I,  J ) =  DEXP  < -P  < I, J ) /T  > 

10  CONTINUE 

DO  20  1=1, NC 
DO  20  K  =  1, NC 
S(K, I > =0. DO 
DO  20  M=1,NG 

S  <  K,  I  >  =S  <  K,  I )  +QT  ( M,  I )  *TAU  ( M,  K  > 
20  CONTINUE 

DO  30  1=1, NC 
F( I ) *1. DO 
DO  30  J=1,NG 

F( I ) *F( I > *QT ( J,  I)*DLOG(S( J, I) ) 
30  CONTINUE 

RETURN 

END 


UVLE  CONTAINS  BUILT-IN  UNIFAC  VLE-PARAMETERS  * 

* 


SUBROUTINE  UVLE ( RI,  QI, AI, WAIN ) 


THE  MAIN  GROUPS  ARE: 


1 

CH2 . 

.  2 

3 

ACH . 

4 

ACCH2  . . . 

.  5 

QH . 

6 

CH30H  .... 

.  7 

H20 . 

8 

ACOH  .... 

9 

CH2C0  . . . 

10 

CHO . 

11 

CCOO . 

12 

HCOQ  .... 

13 

CH20  .... 

14 

CNH2  .... 

15 

CNH . 

16 

<  C ) 3N  .... 

17 

ACNH2  . . . 

18 

PYRIDINE 

19 

CON  • • • • • 

20 

COOH  .... 

21 

CCL . 

22 

CCL2  m m m m 

23 

CCL3  .... 

24 

CCL4  .... 

25 

ACCL  • • • • 

26 

CN02 . 

27 

ACN02  . . . 

28 

CS2 . 

29 

CH3SH  . . . 

30 

FURFURAL 

31 

DOH . 

32 

I  . 

33 

BR . 

34 

C=-C  .... 

35 

DMSO  .... 

36 

ACRY . 

37 

CLCC  .... 

38 

ACF . 

39 

DWF . 

40 

CF2 . 

THE  SUB  GROUPS  ARE: 


1 

CH3 . 

2 

CH2 . 

.  3 

CH . 

4 

c . 

5 

CH2»CH 

6 

CH=CH  .... 

7 

CH2=C  . . . 

.  a 

CH*C  (if* 

.  9 

C*C  • • • • • 

10 

ACH  . . . 

11 

AC  ••••••« 

12 

ACCH3  . . . 

13 

ACCH2  . . . 

14 

ACCH  .... 

15 

OH  . . . . 

16 

CH30H  .... 

17 

H20 . 

18 

ACOH  .... 

19 

CH3C0  . . . 

20 

CH2C0  . 

21 

CHO . 

22 

CH3C00  , . 

23 

CH2C0Q  . . 

24 

HCOO  .... 

25 

CH30  . . 

26 

CH20 . 

27 

CH-Q  .  .  .  . 

28 

FCH20  . . . 

29 

CH3NH2  . . 

30 

CH2NH2 

31 

CHNH2  .... 

32 

CH3NH  . . . 

33 

CH2NH  . . . 

34 

CHNH  .... 

35 

CHON  . . 

36 

CH2N . 

37 

ACNH2  . . . 

38 

C5H5N  . . . 

39 

C5H4N  . . . 

40 

C5H3N  . 

41 

CH3CN  .... 

42 

CH2CN  . . . 

43 

COOH  .... 

44 

HCOOH  . . . 

45 

CH2CL  . 

46 

CHCL . 

47 

CCL  • • • • • 

48 

CH2CL2  . . 

49 

CHCL2  . . . 

50 

CCL2  . . 
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c 

51 

CHCL3  .... 

52 

c 

56 

CH2N02  . . . 

57 

c 

61 

CH2SH  .... 

62 

c 

66 

CH=-C  .... 

67 

c 

71 

ACF . 

72 

c 

76 

CP  ... _ 

c 

c 

c 


CCL3  .... 

53 

CCL4  .... 

54 

CHN02  . . . 

58 

ACN02  . . . 

59 

FURFURAL 

63 

( CH20H ) 2 

64 

C=-C _ 

68 

DMSO  .... 

69 

DMF-1  ... 

73 

DMF-2  . . . 

74 

ACCL  *  • *  * 

55 

CH3N02  . 

CS2  ..... 

60 

CH3SH  . . 

I . 

65 

BR . 

ACHY  .... 

70 

CL  <  C= -C ) 

CF3 . 

75 

CF2  . . . . 

c 


IMPLICIT  REAL*S(A-H, 0-Z) 

COMMON/UNIF/RT  (10#  10),QT(10,  10),TAU<10,  10),S<10,  10  ) ,  F  ( 10  > ,  G  <  10  > , 
R< 10), P< 10, 10) 

DIMENSION  AI <  40,  40 ) ,  RR  <  7S ) ,  GQ ( 76 ) , RI <  76 ) ,  QI ( 76 ) ,  MAINSG ( 76 ) , 

MAIN (76) 


DATA  MAINSG/  4*1,  5*2,  2*3,  3*4,  1*5,  1*6,  1*7,  1*8,  2*9,1*10, 

*  2*11, 1*12, 4*13, 3*14, 3*15,  2*16,  1*17,  3*18,  2*19,  2*20, 

*  3*21, 3*22, 2*23, 1*24, 1*25, 3*26,  1*27,  1*28,  2*29, 1*30, 

*  1*31, 1*32, 1*33, 2*34, 1*35, 1*36, 1*37, 1*38, 2*39, 3*40/ 

C 

DATA  RR/O.  9011,  0.  6744,  O.  4469,  0.  2195,  1.  3454,  1.  1167,  1.  1173,  0.  8886, 
*0.  6605,  0.  5313,  O.  3652,  1.  2663,  1.  0396,  0.  8121,  1.  0,  1.  4311,  O.  92,  O.  8952, 
*1.  6724,  1.  4457,  O.  9980,  1,  9031,  1.  6764,  1. 2420,  1.  1450,  0.  9183,  0.  6908, 
*0.  9183,  1.  5959,  1.  3692,  1.  1417,  1.  4337,  1. 2070,  0.  9795,  1.  1865,  O.  9597, 
*1. 0600,  2.  9993,  2.  8332,  2.  6670,  1.  8701,  1.  6434,  1.  3013,  1.  5280,  1.  4654, 
*1.  2380,  1.  0060,  2.  2564,  2.  0606,  1,  8016,  2.  8700,  2.  6401,  3.  3900,  1.  1562, 
*2.  0086,  1.  7818,  1.  5544,  1.  4199,  2.  0570,  1. 8770,  1.  6510,  3.  1680,  2.  4088, 
*1.  2640,  0.  9492,  1.  2920,  1.  0613,  2.  8266,  2.  3144,  0.  7910,  0.  6948,  3.  0856, 
*2. 6322,  1.  4060,  1.  0105, 0. 6150/ 

C 


DATA  QQ/O.  848,  0.  540,  0.  228,  0.  000,  1.  176,  O.  867,  0.  9 88,  0.  676,  O.  485, 
*0.  400,  0.  120,  0.  968,  0.  660,  0.  348,  1.  200,  1.  432,  1.  400,  0.  680,  1.  488, 
*1.  180,  0.  948,  1.  728,  1.  420,  1.  188,  1.  088,  0.  780,  0.  468,  1.  100,  1.  544, 
*1.  236,  0.  924,  1.  244,  0.  936,  0.  624,  0.  940,  0.  632,  0.  816,  2.  113,  1.  833, 
*1.  553,  1.  724,  1.  416,  1.  224,  1.  532,  1.  264,  0.  952,  0.  724,  1.  988,  1.  684, 
*1,  448,  2.  410,  2.  184,  2.  910,  0.  844,  1.868,  1.  560,  1.  248,  1.  104,  1.  650, 
*1.  676,  1.  368,  2.  481,  2.  248,  0.  992,  0.  832,  1.  088,  0.  784,  2.  472,  2.  052, 
*0.  724,  0.  524,  2.  736,  2.  120,  1.  380,  0.  920,  0.  460/ 


OPEN ( 9, FILE* ' AVLE. DAT ' ,  STATUS* ' OLD ' ,  ACCESS* ' SEQUENTIAL ' , 
♦FORM* *  FORMATTED ' ) 

OPEN ( 10, FILE* ’ AVLE.  OUT ' ,  STATUS* ' NEW ' ,  ACCESS* ' SEQUENTIAL ' , 
♦FORM* ' FORMATTED ' ) 

C 

DO  5  1*1,76 
RI  ( I ) *0. DO 
QI < I ) *0. DO 
MAIN  < I ) *0. DO 
£  CONTINUE 
C 

DO  20  1*1, 40 
DO  10  J=l, 40 

AI ( I,  J) *0. DO 
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nnnnnnnnn  n  non 


10  CONTINUE 

20  CONTINUE 
C 

DO  40  1  =  1,40 
DO  30  J=l, a 

READ  (9,  *,  END=40)  (AI<  ( S* ( J-i >  +K  > ,  I),  K=1,S> 
30  CONTINUE 

40  CONTINUE 
C 

DO  50  1  =  1,76 
RI < I ) =  RR  < I ) 

QI < I ) =QQ( I ) 

MAIN  < I ) =MAINSG < I ) 

50  CONTINUE 
C 

DO  70  1=1, 40 
DO  60  J=l,8 

WRITEdQ,  999)  (  AI  <  (  5*  ( J -1  )  +  K  ) ,  I ) ,  K  =  l,5) 

60  CONTINUE 
70  CONTINUE 

-  FORMAT  STATEMENTS  - 

999  FORMAT (5F10. 3) 

RETURN 
END 


ULLE  CONTAINS  BUILT-IN  UNIFAC  LLE-PARAMETERS 


SUBROUTINE  ULLEt RI,  QI, AI, MAIN ) 

c  »••*•**«»•«•••••*»•••••»••»••••»• 

C 

C  THE  MAIN  GROUPS  ARE : 


C 

1  CH2  . 

C=C . 

3 

ach . 

.  4 

ACCH2  . . . 

.  5 

OH 

C 

6  PI  .  . 

P2 . 

a 

H20 . 

.  9 

ACOH  .... 

10 

L#H3C0  r  1  *  « 

c 

11  CHO  . 

.  12 

FURFURAL 

13 

COOH  .... 

14 

caoc  .... 

15 

CH20 . 

c 

16  CCL  . 

CCL2  *  * • « 

la 

CCL3  • • • • 

19 

CCL4  .... 

20 

ACCL . 

c 

21  CCN  . 

ACNH2  . . . 

23 

CN02  .... 

24 

ACNC2  . . . 

25 

DOH . 

c 

26  DEOH 

PYRIDINE 

28 

TCE . 

29 

MFA . 

30 

DMFA . 

c 

c 

31  TMS  . 

.  32 

DMSO  .... 

c 

c 

THE  SUB 

GROUPS  ARE: 

c 

1  CH3  . 

CH2 . 

3 

CH  . 

.  4 

C . 

.  5 

CH2=CH  . . . 

c 

6  CH=CH 

.  7 

CH  =  C _ 

a 

CH2=C  . . . 

.  9 

ACH . 

10 

AC . 

c 

11  ACCH3 

....  12 

ACCH2  . . . 

13 

ACCH  .... 

14 

OH . 

15 

PI . 
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c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 


c 


c 


c 


c 


16  P2 .  17  H20 .  ia  ACOH _  19  CH3C0  ...  20  CH2C0  .  .  . 

21  CHO  ......  22  FURFURAL  23  COOH  ....  24  HCOOH  ...  25  CH3C00  . . 

26  CH2COO  ...  27  CH30 _  2a  CH20 _  29  CH-0 _  30  FCH20  .  .  . 

31  CH2CL _  32  CHCL _  33  CCL .  34  CH2CL2  .  .  35  CHCL2  .  .  . 

36  CCL2  .  37  CHCL3  ...  38  CCL3  ....  39  CCL4  ....  40  ACCL  _ 

41  CH 3CN  _  42  CH2CN  ...  43  ACNH2  ...  44  CH3N02  , .  45  CH2N02  . . 

46  CHN02  _  47  ACN02  ...  48  (CH20H52  49  (H0MMJ20  50  C5H5N  ... 

51  C5H4N  _  52  C5H3N  ...  53  CCL2=CHCL  54  HC0NHCH3  55  HC0N(CH3> 

56  ( CH214S02  57  (CH2>2S0 


IMPLICIT  REAL*8(A-H, 0-Z) 

C0MM0N/UNIF/RT<10,  10),  QT< 10,  10>,TAU<10,  10 ) ,  S ( 10, 10 > , F < 10 ) , Q < 10 > , 

*  R<10>, P<10, 10) 

DIMENSION  AI ( 40, 40 ) ,  RR < 57 ) ,  QQ <  57 ) ,  RI ( 76 ) ,  01 ( 76 ) , MAINSG  <  57 ) , 

*  MAIN (76) 

DATA  MAINSG/  4*1,  4*2,  2*3,  3*4,  1*5,  1*6,  1*7,  1*8,  1*9,2*10, 

*  1*11,  1*12,  2*13,  2*14,  4*15,  3*16,  3*17,  2*18,  1*19,  1*20, 

*  2*21, 1*22, 3*23, 1*24,  1*25,  1*26,  3*27,  1*28,  1*29,  1*30, 

*  1*31,1*32/ 

DATA  RR/O.  9011,  0.  6744,  O.  4469,  0.  2195,  1.  3454,  1.  1167,  0.  8886,  1.  1173, 
*0.  5313,  0.  3652,  1.  2663,  1.  0396,  0.  8121,  1.  0000,  3.  2499,  3.  2491,  O.  9200, 
*0.  6952,  1.  6724,  1.  4457,  0.  9980,  3.  1680,  1.  3013,  1.  5280,  1.  9031,  1.  6764, 
*1.  1450,  0.  9183,  0.  6908,  0.  9183,  1.  4654,  1.  2380,  1.  0060,  2.  2564,  2.  0606, 
*1.  8016,  2.  8700,  2.  6401,  3.  3900,  1.  1562,  1.  8701,  1.  6434,  1.  0600,  2.  0086, 
*1.  7818,  1.  5544,  I.  4199,  2.  4088,  4.  0013,  2.  9993,  2.  8332,  2.  6670,  3.  3092, 
*2.  4317,  3.  0856,  4.  0358,  2.-8266/ 

DATA  QQ/O.  848,  O.  540,  0.  228,  0.  000,  1.  176,  0.  867,  0.  676,  0.  988,  0.  400, 
*0.  120,  0.  968,  0.  660,  0.  348,  1.  200,  3.  128,  3-.  124,  1.  400,  0.  680,  1.  488, 

*1.  180,  0.  948,  2.  484,  1.  224,  1.  532,  1.  728,  1.  420,  1.  088,  0.  780,  0.  468, 

*1.  100,  1.264,  0.  952,  0.  724,  1.  988,  1.  684,  1.  448,  2.  410,  2.  184,  2.  910, 

*0.  844,  1. 724,  1.  416,  0.  816,  1.  868,  1.  560,  1.  248,  1.  104,  2.  248,  3.  568, 

*2.  113,  1.  833,  1.  553,  2.  860,  2.  192,  2.  736,  3.  200,  2.  472/ 

OPEN (11, FILE= ' ALLE.  DAT ' ,  STATUS® ' OLD ' ,  ACCESS® ' SEQUENTIAL ' , 

*  FORM® ' FORMATTED  * ) 

OPEN ( 12, FILE® ' ALLE. OUT ' ,  STATUS® ' NEW ' , ACCESS® ' SEQUENTIAL ' , 

♦FORM® ' FORMATTED ' > 

DO  5  1=1, 76 
RI ( I ) =0. DO 
QI < I >  =0. DO 
MAIN ( I ) =0. DO 
5  CONTINUE 

DO  20  1=1,40 

DO  10  J= 1, 40 
AI  <  I,  J)  =0.  DO 
10  CONTINUE 


non 


20  CONTINUE 
C 

DO  40  1  =  1,  32 
DO  30  J=i, 6 

READf  11,  »,  END  =  40  )  (  AI  (  <  5*  <  J-l  > +K  >,  I  >,  K  =  I,5> 

30  CONTINUE 

READf  11,  *,  END  =  40)  <AI(K,  I),  K  =  31,32> 

40  CONTINUE 
C 

DO  50  1=1, 57 
RI ( I ) =RR ( I > 

QI ( I > =  00 < I ) 

MAIN ( I ) =MAINSG ( I > 

«  SO  CONTINUE 

C 

DO  70  1=1, 32 
DO  60  J  =  l,  6 

WRITE<12,  999)  (  AI  <  <  5*  (  J-l  )  ♦K  ) ,  I  > ,  K=l,5> 

60  CONTINUE 

WRITE(12,  999)  (AI(K, I),  K  =  3i,32) 

70  CONTINUE 

-  FORMAT  STATEMENTS  - 


C 

C  THE  SUB  GROUPS  ARE: 

C  1  CH3 .  2  CH2 .  3  CH .  4  C .  5  CH2=CH  .  .  . 
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C 

6 

CH*CH 

*  •  •  * 

.  7 

CH2*C  . . . 

.  a 

CH*C  .... 

9 

C  =  C . 

10 

ACH  .... 

c 

11 

AC  .  .  . 

•  •  *  « 

12 

ACCH3  . . . 

13 

ACCH2  . . . 

14 

A CCH  .... 

15 

OH . 

C 

16 

CH30H 

•  a  a  a 

17 

H20 . 

18 

ACOH  .... 

19 

CH3C0  . . . 

20 

CH2C0  . . 

C 

21 

CHO  .  . 

*  a  *  * 

22 

CH3C00  . . 

23 

CH2C0Q  . . 

24 

HCOO  ..... 

25 

CH30  . . . 

c 

26 

CH20  . 

•  •  •  a 

27 

CH-0 _ 

28 

FCH20  . , . 

29 

CH3NH2  . . 

30 

CH2NH2  . 

c 

31 

CHNH2 

*  *  *  a 

32 

CH3NH  . . . 

33 

CH2NH  . . . 

34 

CHNH  .... 

35 

CH3N  . . . 

c 

36 

CH2N  . 

•  *  *  * 

37 

ACNH2  . . . 

38 

C5H5N  . . . 

39 

C5H4N  . . . 

40 

C5H3N  . . 

c 

41 

CH3CN 

•  •  «  • 

42 

CH2CN  . . . 

43 

COOH  .... 

44 

HCOOH  . . . 

45 

CH2CL  . . 

c 

46 

CHCL  . 

•  •  *  • 

47 

CCL . 

48 

CH2CL2  . . 

49 

CHCL2  . . . 

50 

CCL2  . . . 

c 

51 

CHCL3 

•  •  •  « 

52 

CCL3  . . . 

53 

CCL4 _ 

54 

ACCL  .... 

55 

CH3N02  . 

c 

56 

CH2N02 

57 

CHN02  . . . 

58 

ACN02  . . . 

59 

CS2  ..... 

60 

CH3SH  . . 

c 

61 

CH2SH 

•  •  •  • 

62 

FURFURAL 

63 

( CH20H ) 2 

64 

X  ••••••• 

65 

BR . 

c 

66 

CH= -C 

•  «  •  • 

67 

C=~C 

68 

DMSO  .... 

69 

ACRY  .... 

70 

CL<C*-C> 

c 

71 

ACF  .  . 

•  *  «  « 

72 

DMF-i  ... 

73 

DMF-2  ... 

74 

CF3 . 

75 

CF2  .... 

C  76  CF 

C 

. . . 

c 


c 


c 


IMPLICIT  REAL*8< A-H, 0-2) 

COMMON/ UNIF/RT( 10,  10) , QTC 10, 10) , TAUC 10,  10 ) ,  S < 10,  10 > ,  F ( 10 > , Q < 10 > , 
R( 10),  P< 10, 10) 

DIMENSION  AI <  40, 40 ) , RR  <  76  > , GO <  76 ) , RI <  76  > , QI <  76 ) , MAINSG  <  76  > , 

MAIN (76) 

DATA  MAINSG/  4*1,  5*2,  2*3,  3*4,  1*5,  1*6,  1*7,  1-8,  2*9,1*10, 
2*11,  1*12,  4*13,  3*14,  3*15,  2*16, 1*17,  3*18,  2*19,  2*20, 
3*21,  3*22,  2*23,  1*24,  1*25,  3*26,  1*27,  1*28,  2*29,  1*30, 
1*31,  1*32,  1*33,2*34,  1*35,  1*36,  1*37,  1*38, 2*39, 3*40/ 


DATA  RR/O.  9011,  0.  6744,  0.  4469,  0.  2195,  1.  3454,  1,  1167,  1.  1173,  0.  8886, 
*0.  6605,  0.  5313,  0.  3652,  1.  2663,  1 . 0396,  0.  8121,  1.  0,  1. 4311,  0.  92,  0.  8952, 
*1.  6724,  1.  4457,  O.  9980,  1. 9031,  1.6764,  1.  2420,  1.  1450,  0.  9183,  0.  6908, 
•0.  9183,  1.  5959,  1.  3692,  1.  1417,  1.  4337,  1.  2070,  0.  9795,  1.  1865,  0.  9597, 

*1.  0600,  2.  9993,  2.  8332,  2.  6670,  1. 8701,  1.  6434,  1.  3013,  1.  5280,  1.  4654, 

*1.  2380,  1.  0060,  2.  2564,  2.  0606,  1.8016,  2.  8700,  2.  6401,  3.  3900,  1.  1562, 
*2.  0086,  1.  7818,  1.  5544,  1.  4199,  2.  0570,  1. 8770,  1.  6510,  3.  1680,  2.  4088, 

*1.  2640,  0.  9492,  1.  2920,  1.  0613,  2.  8266,  2.  3144,  0.  7910,  0.  6948,  3.  0856, 

*2. 6322,  1. 4060,  1. 0105,  0. 6150/ 


DATA  QQ/Q.  848,  0.  540,  O.  228,  0.  000,  1 .  176,  0.  867,  0.  9 88,  0.  676,  O.  485, 
*0.  400,  0.  120,  0.  968,  0.  660,  0.  348,  1 . 200,  1.  432,  1.  400,  0.  680,  1.  488, 

*1.  180,  0.  948,  1.  728,  1.  420,  1.  188,  1 . 088,  0.  780,  0.  468,  1.  100,  1.  544, 

*1.  236,  0.  924,  1. 244,  0.  936,  0.  624,  0.  940,  0.  632,  0.  816,  2.  113,  1.  833, 
*1.  553,  1.  724,  1. 416,  1.  224,  1.  532,  1. 264,  O.  952,  0.  724,  1.  988,  1.  684, 
*1.  448,  2.  410,  2.  184,  2.  910,  0.  844,  1. 868,  1.  560,  1.  248,  1.  104,  1.  650, 

*1.  676,  1.  368,  2.  481,  2.  248,  0.  992,  0.  832,  1.  088,  0.  784,  2.  472,  2.  052, 

•0.  724,  0.  524,  2.  736,  2.  120,  1.  380,  0.  920,  O.  460/ 


OPEN ( 13, FILE= 'AENV. DAT' , STATUS* 'OLD',  ACCESS* 'SEQUENTIAL', 
-FORM* 'FORMATTED' ) 

OPEN ( 14, FILE* 'AENV. OUT ' , STATUS* ' NEW ' ,  ACCESS* 'SEQUENTIAL', 
-FORM* ' FORMATTED ' ) 
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o  n  o 


DO  S  1*1, 76 
RI<  I  > *0. DO 
OKI)  *0.  DO 
MAIN ( I ) *0. DO 
5  CONTINUE 

DO  20  1  =  1,  40 
DO  10  J*i, 40 

AKI,  J)  =0.  DO 
10  CONTINUE 
20  CONTINUE 

DO  40  1*1,40 
do  30  j*:,a 

READ <13,  »,  END=40>  <  Al  <  <  5*  <  J-i  )  «-K >  ,  I  >  ,  K*i,5> 
30  CONTINUE 

40  CONTINUE 

DO  50  1  =  1,76 
RI ( I ) =RR  < I ) 

QICI) =Q0< I ) 

MAIN  < I ) =MAINSG< I ) 

50  CONTINUE 

DO  70  1*1,40 
DO  60  J*l,8 

WRITE  <  14,  999  >  <AI<  <5*<J-i)*K),  I),  K*l,5) 

60  CONTINUE 

70  CONTINUE 

-  FORMAT  STATEMENTS  - 

999  FORMAT (5F10. 3) 

RETURN 

END 
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(The  reverse  of  fhis  page  is  blank.) 


